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HARMONIC  EMISSION  FROM  HIGH  POWER 
GYROTRON  OSCILLATORS 

by 

JOSEPH  LAWRENCE  BYERLY 

Submitted  to  the  Department  of  Physics 
on  May  11,  1984  in  partial  fulfillment  of  the  requirements 
for  the  Degree  of  Master  of  Science  in  Physics 

ABSTRACT 

V  Experimental  and  theoretical  results  are  presented  for  the  second  harmonic  emission 
(-•  %  2^7.  where  u>c  is  the  electron  cyclotrpnJreguencyffrom  a  pulsed,  high  power  gyrotron 
designed  for  high  frequency  fundamental  («;  ss  1*7)  operation.  The  gyrotron  produces  more 
than  100  kW  of  output  power  at  140  GHz  in  the  TEo.3,1  mode  at  the  fundamental.  Eight 
second  harmonic  modes  were  observed  in  gyrotron  operation  between  4.0  and  6.0  tesla 
with  65  kV  beam  voltage  and  up  to  five  amps  of  beam  current.  The  two  strongest  modes 
were  the  TEn.2.1  at  241  GHz  which  produced  approximately  25  kW  of  output  power  at 
7('c  efficiency  and  the  TE92.1  with  ss  15  kW  at  209  GHz.  These  represent  the  highest 
power  measured  for  high  frequency  (>  100  GHz)  gyrotron  operation  at  a  harmonic. 

The  theoretical  investigation  involved  a  linear  theory  for  determining  threshold  oscil¬ 
lation  conditions  and  two  nonlinear  theories  for  efficiency  calculations.  For  the  gyrotron 
resonator  which  is  excited  at  the  fundamental  in  the  TEo  3,1  mode  at  140  GHz,  there  are 
36  different  modes  in  the  frequency  range  203  303  GHz  which  could  be  excited  atu  «  2 u>c 
for  a  magnetic  field  between  4.0  6.0  tesla.  The  linear  theory  results  predicted  28  of  these 
second  harmonic  resonator  modes  to  be  above  threshold  for  operation  at  65  kV  and  five 
amps.  Minimum  starting  currents  for  the  strongest  second  harmonic  modes  were  calcu¬ 
lated  to  be  0.5  A.  which  is  of  the  same  order  as  starting  currents  for  fundamental  modes. 
The  calculations  using  both  slow  and  fast  time  scale  nonlinear  theories  predicted  overall 
efficiencies  of  ss  20°A.  This  is  only  a  factor  of  two  lower  than  the  theoretical  fundamen¬ 
tal  mode  efficiencies.  The  relatively  high  second  harmonic  efficiencies  arise  because  these 
modes  have  higher  Q  factors  which  compensate  for  the  reduced  gain  at  the  harmonics. 

The  measured  harmonic  emission  was  less  than  the  theoretical  predictions  both  in 
terms  of  number  of  modes  observed  and  the  overall  efficiencies.  A  number  of  possible 
explanations  for  the  discrepancy  were  considered.  At  most  magnetic  field  values  there 
was  some  radiation  at  the  fundamental  and  the  second  harmonic  signal  was  weak  or 
nonexistent.  The  two  strong  second  harmonic  modes  appeared  at  magnetic  field  values 
where  there  was  no  fundamental  oscillation.  This  seemed  to  indicate  that  oscillation  at 
the  fundamental  was  suppressing  the  weaker  and  missing  second  harmonic  modes.  A 
qualitative  consideration  of  the  mode  competition  between  the  fundamental  and  second 
harmonic  supports  this  interpretation  of  mode  suppression. 

The  present  study  indicates  that  second  harmonic  emission  can  achieve  high  power 
levels  at  very  high  frequencies  with  reasonable  efficiencies.  Using  available  superconducting 
magnets  and  second  harmonic  operation,  high  power  emission  well  into  the  submillimeter 
(  >  600  GHz)  appears  feasible 
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Chapter  1 

Introduction 

Within  the  past  two  decades,  the  gyrolron  has  rapidly  developed  as  an  efficient 
source  of  high  frequency  microwave  radiation.  A  worldwide  effort  of  experimental 
gyrotron  designs  has  produced  a  wide  range  of  operating  parameters  including 
output  power  levels  in  excess  of  a  megawatt,  pulse  lengths  from  microseconds 
to  continuous  (cw),  and  wavelengths  from  several  centimeters  to  well  into  the 
submillimeter.  The  technology  and  demand  for  gvrotrons  is  such  that  they  are 
now  available  commercially. 

Increasing  with  the  performance  level  of  gvrotrons  is  the  list  of  proposed 
applications.  Among  these  are  high  frequency  radar,  particle  accelerators,  spec¬ 
troscopy  of  biological  materials,  and  as  a  scattering  source  for  plasma  diagnos¬ 
tics.  But  the  most  impressive  application  thus  far  has  been  in  electron  cyclotron 
resonance  heating  (ECRH)  experiments  in  magnetically  confined  fusion  devices. 
Significant  bulk  plasma  heating  increases  have  been  recorded  in  the  Soviet  Union 
1  and  gvrotrons  for  ECRU  have  become  an  integral  part  of  many  experimental 
fusion  machines  in  the  United  States.  The  ECRH  demands  of  a  commercially 


viable  fusion  reactor  have  driven  the  long  term  gyrotron  goals  to  long  pulse,  high 
efficiency,  with  megawatt  power  levels  in  the  100-300  GHz  frequency  range.  As 
stated  previously,  each  of  these  goals  alone  has  been  realized,  but  no  gyrotron 
exists  today  which  incorporates  all  of  them. 

1.1  —  Historical  Development  of  the  Gyrotron 

Microwave  devices  in  the  1950’s  all  had  similar  limitations  which  precluded 
design  and  operation  in  the  millimeter  and  submillimeter  regime.  Internal  com¬ 
ponent  dimensions  were  directly  related  to  the  wavelength  (A)  of  the  output  radi¬ 
ation.  Systems  such  as  the  travelling  wave  tube  (TYYT)  or  magnetron  required  a 
periodic  structure  on  the  order  of  A  and  klystrons  had  electrodes  with  dimensions 
less  than  A.  Obviously,  radiation  wavelengths  of  a  few  millimeters  or  less  required 
a  new  mechanism  for  amplification. 

Such  a  mechanism  was  first  theorized  in  1958  by  R.  Q.  Twiss  2  .  An  astro¬ 
physicist.  Twiss  surmised  that  under  certain  conditions,  electromagnetic  radiation 
traversing  a  magnetized,  relativistic  plasma  would  be  amplified.  In  1959.  Schneider 
3  using  (piant urn  mechanics  and  A.  V.  Gaponov  •)  employing  a  classical  treat¬ 
ment.  separately  reached  a  similar  conclusion  concerning  the  possibility  of  gain 
from  relativistic  electrons  gyrating  in  a  magnetic  field.  H irsh field  and  Wachlel  5 
provided  the  first  experimental  confirmation  of  the  gain  mechanism  in  1961. 

The  amplified  radiation  from  the  theorized  interaction  was  necessarily  near 
the  Doppler  shifted  cyclotron  resonance.  With  ^  representing  the  angular  fre¬ 
quency  and  k  the  wave  vector  of  the  radiation,  the  condition  may  be  stated 


where  n  is  an  integer  representing  the  harmonic  number,  v  the  drift  velocity  and 
uc  is  the  relativistic  electron  gyrofrequency.  That  is 


eB 

u )  c  — 

irne 


(1.2) 


with  e  and  me  the  electron  charge  and  rest  mass,  B  the  magnetic  field  strength 
and  7  is  the  relativistic  factor  given  by 


1 


(1.3) 


where  J  -  r  c  is  the  velocity  normalized  to  the  speed  of  light  (c).  Appropri¬ 
ately  enough,  devices  employing  this  amplification  mechanism  were  called  electron 
cyclotron  masers  (ECM)  or  cyclotron  resonance  masers  (CRM),  the  latter  term 
appearing  extensively  in  Soviet  literature. 

Theoretical  and  experimental  work  continued  throughout  the  1960’s  with  a 
massive  effort  taking  place  in  the  Soviet  Union.  Many  possible  configurations  of 
electron  beam  and  magnetic  field  were  considered.  However,  the  most  promising 
was  a  helical,  annular  electron  beam  from  a  magnetron  injection  gun  transiting  an 
open  resonator  near  cutoff,  aligned  with  a  static  magnetic  field.  This  combined 
t  he  capability  to  confine  an  intense,  monocnergetic  electron  beam  with  the  decided 
advantage  of  simplicity  of  design  and  construction.  By  operating  near  the  cutoff 
frequency  of  the  resonator  cavity,  the  k  -v  term  is  small  so  that  high  efficiency  can 
be  achieved  with  electron  beams  having  a  large  velocity  spread.  The  term  gvrotron 
was  first  used  by  Gaponov  [6  for  devices  of  this  description.  The  name  gvrotron 
has  since  acquired  the  additional  restrictive  meaning  of  weakly  relativistic  electron 
beam  voltages,  such  that  02  <:  1. 


9 


In  the  latter  half  of  the  1960’s,  interest  in  this  class  of  free  electron  device 
was  minimal  in  the  west,  but  the  Soviet  effort  continued  with  significant  theo¬ 
retical  advances.  In  particular,  a  nonlinear  theory  developed  Gaponov,  Petelin, 
and  Yulpatov  |7j  and  reduced  to  a  particularly  tractable  form  by  Nusinovich  and 
Erm  [8] ,  allowed  gyrotron  efficiency  calculations  and  thus  optimization  of  device 
parameters  including  resonator  cavity  profile.  The  results  of  the  theoretical  work 
became  evident  in  the  early  1970’s  when  Soviet  experimentalists  provided  impres¬ 
sive  milestones  of  achievement  .  These  results  sparked  renewed  interest  in  gyrotron 
research  in  the  west. 

A  diagram  of  the  status  of  gyrotron  research  to  date  is  shown  in  Figures  1.1a 
and  1.1b.  Figure  1.1a  considers  only  gvrotrons  operating  at  the  fundamental 
(n  ~  1)  while  Figure  1.1b  depicts  the  results  at  higher  harmonics  (n  >  1).  Both 
graphs  show  achievements  in  terms  of  frequency  and  output  power  with  separate 
symbols  for  various  degrees  of  pulsed  operation.  For  these  figures,  operation  with 
pulse  duration  of  less  than  10  msec  was  considered  short  pulse,  while  greater  than 
50  msec  is  shown  as  long  pulse. 

Examination  of  Figure  1.1a  reveals  the  tradeoff  between  long  pulse  length  and 
high  output  power  at  high  frequency.  The  highest  outpul  power  in  cw  operation 
is  the  3-12  kW  Varian  [9j  result  at  28  GHz.  Varian  j  1 0j  and  Hughes  11  have 
both  produced  200  kW,  60  GHz  gvrotrons.  The  Hughes  version  is  long  pulsed  and 
the  Varian  design  is  capable  of  cw  operation.  It  should  be  mentioned  that  chart 
space  in  Figure  1.1a  did  not  allow  inclusion  of  the  increasing  number  of  commercial 
ventures  in  this  low  frequency,  long  pulse  regime.  In  particular,  Thomson-CSF  [12; 
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and  Nippon  Electric  Company  (NEC)  [13]  have  produced  ~  200  kW  gyrotrons  for 
long  pulse  operation  at  35  GHz.  Toshiba  has  also  constructed  22  and  70  GHz 
gyrotrons  [14].  The  Soviet  contribution  in  this  area  of  operation  are  gyrotrons 
presently  used  to  heat  the  Kurchatov  T-10  Tokamak  [15].  These  devices  operate 
at  85  GHz  with  an  average  output  power  of  200  kW  in  pulses  of  up  to  about  50 
msec. 

Increasing  the  frequency  of  cw  operation  has  resulted  in  significant  decreases 
in  the  output  power  available.  As  shown  in  Figure  1.1a.  the  highest  frequency 
yet  achieved  in  cw  operation  at  the  fundmental  is  150  GHz  with  22  k\V  of  output 
power  reported  by  Andronov  et  al.  16  .  This  surpassed  the  1974  result  described 
by  Zaytsev  17  of  a  cw  gyrotron  capable  of  12  kW  at  107  GHz. 

Short  pulsed  operation  greatly  increases  the  operating  space,  both  in  terms 
of  power  and  frequency.  The  1  MW  power  level  was  first  exceeded  by  a  gyrotron 
operating  in  ^hort  pulses  at  the  fundamental  in  1978  as  reported  by  Andronov 
1G  .  These  results  included  the  achievement  of  1.25  MW'  at  45  GHz  and  1.1  MW 
at  100  GHz.  Experiments  by  GoFdenberg  which  were  reported  by  Flyagin  18 
recent  l\  exceeded  this  power  level.  Using  a  90  kY.  50  A  eleclron  beam  with  a 
10  psec  pulse  length.  1.4  MW'  of  output  power  was  achieved  at  100  GHz.  The 
same  paper  described  experiments  well  into  the  submillimeter  by  Luchinin  and 
Nusinovich.  which  are  represented  in  Figure  1.1a  by  the  six  points  beyond  300 
GHz  with  power  levels  from  20  to  130  kW.  These  dramatic  results  were  produced 
by  using  pulsed  magnets  capable  of  23  tesla. 
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The  final  points  in  Figure  1.1a  represent  short  pulsed  results,  but  were  de¬ 
signed  for  extension  to  the  cw  and  ECRH  applications.  Wall  loading  effects  on  the 
resonator  cavity  were  considered  and  continuous  magnets  were  used.  A  gyrotron 
at  the  Naval  Research  Laboratory  (NRL)  [19]  achieved  340  kW  of  output  power 
at  35  GHz  with  an  incredible  53%  efficiency.  The  MIT  gyrotron  [20]  produced  175 
kW  of  power  at  140  GHz  in  operation  at  the  fundamental.  This  represents  the 
highest  power  at  high  frequency  (>  100  GHz)  in  a  low  ohmic  loss  resonator  mode 
which  is  compatible  with  scaling  up  to  long  pulse  or  cw  operation. 

High  power  gyrotrons  operating  at  harmonics  of  the  cyclotron  frequency  are 
shown  in  Figure  1.1b.  All  of  the  results  depicted  are  second  harmonic  (n  -  2) 
except  for  one  third  harmonic  result  which  is  labelled. 

A  large  part  of  the  Soviet  experimental  effort  in  the  early  1970’s  used  gy- 
rotrons  designed  for  operation  at  the  second  harmonic.  In  1974,  Rise!  ]2l]  ob¬ 
tained  10  kW  of  cw  power  and  30  kW  in  pulsed  operation  at  34  GHz.  The  next 
month,  Zaytsev  j  1 7 .  published  results  of  7  kW  at  154  GHz  in  pulsed  operation, 
with  2.4  kW  at  157  GHz  and  1.5  kW  at  326  GHz  in  cw  operation.  Recent  Soviel 
harmonic  gyrotron  investigations  by  Zapevalov  reveal  a  marked  increase  in  power 
and  efficiency  by  using  a  stepped  profile  cavity.  In  1979  ,22  .  pulsed  operation 
at  25  GHz  produced  500  kW.  The  latest  Zapevalov  results  (23]  include  a  pulsed 
second  harmonic  gyrotron  emitting  150  kW  of  output  power  at.  34  GHz.  and  third 
harmonic  pulsed  operation  producing  100  kW  of  power  at  55  GHz.  This  last  result 
represents  a  remarkable  10%  efficiency  at  the  third  harmonic. 


Fig  1.1b  —  Gyrotrons  Operating  at  Higher  Harmonics 
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There  have  been  few  investigations  of  high  power,  high  frequency  harmonic 
operation  outside  of  the  Soviet  Union.  The  Harry  Diamond  Laboratories  in  con¬ 
junction  with  NRL  [24'  completed  a  design  effort  for  a  240  GHz,  ~  10  kW  gyrotron, 
but  initial  experimental  results  have  yet  to  produce  the  expected  efficiencies.  Some 
work  has  been  carried  out  in  the  People’s  Republic  of  China  (PRC),  including  30 
kW  of  power  at  37  GHz  [25].  Thomson-CSF  '12  has  recently  reported  on  initial 
studies  of  a  second  harmonic  gyrotron  which  has  produced  ~  30  kW  at  70  GHz. 
The  only  other  high  power,  second  harmonic  operation  of  a  gyrotron  is  the  MIT 
result  : 26] ,  which  is  the  basis  for  this  thesis.  As  depicted  in  Figure  1.1b,  the  25  kW 
of  output  power  at  241  GHz  achieved  in  the  present  study  is  a  singular  achieve¬ 
ment  for  harmonic  operation.  In  the  future,  it  may  be  possible  to  extend  these 
results  to  higher  power  and  frequency. 

1.2  -  The  MIT  Gyrotron 

The  MIT  results  cited  in  the  previous  section  were  produced  by  a  gyrotron 
which  has  been  in  operation  since  1982  al  the  Francis  Bitter  National  Magnet 
Laboratory  on  the  MIT  campus.  The  design  and  intitial  results  were  reported  in 
27  .  Recent  results  describing  enhanced  efficiency  and  improved  diagnostics  are 
found  in  20  . 

Figure  1.2  is  a  schematic  diagram  of  the  device.  The  overall  axisyrnmetric 
geometry  is  typical  of  most  gyrotron  designs.  To  the  left  in  Figure  1.2  is  the 
electron  gun  which  produces  an  annular  electron  beam.  The  beam  travels  to  the 
right  entering  the  magnetic  compression  region  and  through  the  resonator  cavity. 
Exiting  the  resonator  and  high  magnetic  field  region,  the  beam  is  deposited  on 


the  collector  which  also  serves  as  an  oversized  output  waveguide  for  the  generated 
RF.  The  microwave  radiation  passes  through  a  fused  quartz  output  window  which 
allows  the  entire  resonator  tube  section  to  be  kept  under  vacuum,  normally  less 
than  4  x  10~8  torr. 

The  electron  gun  is  a  magnetron  injection  type  constructed  by  Varian  As¬ 
sociates  j28'  which  produces  a  nonlaminar  flow  of  electrons.  The  designed  ratio 
of  perpendicular  to  parallel  velocities  {v± /v\\)  with  respect  to  the  main  magnetic 
Held  is  1.49.  Beam  cathode  voltages  of  up  to  80  kV  have  been  used  producing 
beam  currents  of  up  to  10  A. 

The  Bitter  magnet  which  produces  the  main  solenoidal  field  is  a  copper  core, 
water  cooled  magnet  capable  of  generating  fields  in  excess  of  10  T  in  cw  operation. 
Additionally,  two  coils  capable  of  generating  0.2  T  are  located  symmetrically  about 
the  electron  gun  cathode  region  to  allow  minor  changes  in  the  magnet  ic  compres¬ 
sion  and  beam  quality.  These  are  labelled  ‘'gun  coils"  in  Figure  1.2. 

The  resonator  cavity  is  a  section  of  cylindrical  waveguide  with  a  downlaper 
on  the  beam  input  end  and  an  uptaper  to  the  output  waveguide  on  the  other  end. 
1  he  input  section  downlaper  is  intended  to  keep  the  microwave  radiation  from 
leaking  back  towards  the  electron  gun. 

The  same  resonator  produced  both  the  140  GHz.  175  kW  fundamental  result 
and  the  241  GHz,  25  kW  second  harmonic  result  shown  in  Figures  1.1a  and  1.1b. 
The  cavity  was  designed  for  fundamental  operation  in  the  TE,  1,3.1  mode  with  a 
straight  section  length  of  6A  and  a  radius  of  0.348  cm.  The  theoretical  diffractive  Q 
(Qi>)  of  the  cavity  was  1500  at  the  fundamental.  For  cw  applicability,  the  cavity 
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Fig  1.2  —  Schematic  Diagram  of  the  MIT  Gyrotron 
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wall  loading  was  constrained  to  2  kW /cm2  and  the  electron  beam  was  placed  to 
interact  with  the  second  radial  maximum  of  the  electromagnetic  field.  This  unique 
feature  reduces  space  charge  forces  in  the  beam  which  could  degrade  efficiency. 

The  gyrotron  operates  in  a  pulsed  mode  with  a  pulse  length  of  1  /zsec  and  a 
repetition  rate  of  4  Hz.  Nominal  beam  voltage  and  current  for  the  investigation 
described  by  this  paper  were  64.3  kV  and  5  A.  The  electron  beam  guiding  center 
radius  in  the  resonator  was  designed  for  0.182  cm.  but  adjustments  to  the  gun  coils 
affecting  the  magnetic  compression  yield  0.197  cm  as  a  more  realistic  figure.  These 
values  are  mentioned  now,  since  they  are  important  parameters  in  the  theories 
discussed  in  Chapter  2. 

1.3  —  Motivation  for  an  Harmonic  Emission  Investigation 

Harmonic  emission  from  high  power  gyrotrons  demands  investigation  from 
two  overall  distinct  points  of  view.  The  first  is  that  gyrotrons  designed  for  har¬ 
monic  operation  are  viable  for  some  applications.  The  second  is  that  harmonic, 
emission  can  produce  serious  effects  on  a  gyrotron  intended  for  operation  at  the 
fundamental. 

Examining  these  in  order,  Soviet  experimentalists  have  designed  gyrotrons  for 
harmonic  operation  and  achieved  impressive  results.  Analytic  theories  discussed  in 
Chapter  2  predict  reasonable  efficiencies  (  —  20%)  at  the  second  harmonic  for  a  gy¬ 
rotron  operating  with  the  parameters  of  the  present  device.  The  major  advantage 
of  harmonic  operation  is  the  lowering  of  the  magnet  ic  field  requirement.  Since  in 
a  gyrotron,  uj  nwL  and  ~  D ,  going  to  second  harmonic  operation  effectively 
halves  the  magnetic  field  required  to  produce  the  same  frequency.  Efficient,  high 


harmonic  operation  therefore  leads  to  the  possibility  of  using  permanent  magnets 
at  some  time  in  the  future. 

Achieving  high  efficiency  in  high  harmonic  operation  was  the  goal  of  three 
recent  studies.  Experiments  by  Hirshfield  |29]  on  a  CRM  employing  a  solid  beam 
Pierce  gun  have  produced  emission  at  up  to  the  ninth  harmonic  at  wavelengths 
from  8  mm  to  less  than  1  mm.  The  theoretical  interpretation  presented  by  the 
same  author  in  an  earlier  paper  [30]  and  supported  by  these  results,  is  that  cou¬ 
pling  at  the  beam  boundary  excites  Bernstein  modes  on  the  beam.  Lau  |3]  has 
proposed  a  device  called  a  gvro-magnetron  which  has  a  corrugated  resonator  wall 
similar  to  a  magnetron.  Emission  would  occur  predominantly  at  the  cyclotron 
harmonic  corresponding  to  the  number  of  vanes  in  the  cavity.  A  mildly  relativis¬ 
tic,  annular  electron  beam  would  still  be  employed.  Finally,  McDermott  et  al.  [32] 
have  reported  on  a  gyrotron  with  a  unique  axis-encircling  electron  beam.  Initial 
results  have  produced  multi-kW  power  levels  from  the  eleventh  harmonic  at  65 
GHz. 

The  second  major  reason  for  studying  harmonic  emission  is  to  understand 
its  effects  on  a  gyrotron  designed  for  operation  at  the  fundamental.  I'nrest rained 
harmonic  oscillation  could  possibly  damage  tube  components  such  as  the  output 
window.  Mode  competition  between  the  desired  fundamental  mode  and  parasitic 
harmonic  modes  could  result  in  degraded  efficiency. 

Finally,  megawatt  designed  gyrotrons  in  the  future  will  necessitate  the  use  of 
highly  overmoded  resonators.  A  study  of  the  harmonic  emission  permits  the  inves- 
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tigation  of  an  overmoded  cavity.  For  example,  the  MIT  gyrotron  has  a  resonator 
diameter  of  ~  3A  at  the  fundamental  and  therefore  ~  6A  at  the  second  harmonic. 

This  paper  details  an  harmonic  emission  study  on  the  MIT  gyrotron.  Chap¬ 
ter  2  presents  linear  and  nonlinear  theories  and  applies  them  to  predict  the  second 
harmonic  modes  that  should  be  excited  and  their  efficiency.  Chapter  3  describes 
the  experimental  investigation  to  identify  harmonic  emission,  measure  starting 
current  and  output  power.  Chapter  4  discusses  the  results  and  offers  explanations 
for  discrepancies  between  theory  and  experiment.  Chapter  5  considers  the  pos¬ 
sibility  of  higher  harmonic  emission  with  n  >  2  and  Chapter  6  concludes  with  a 
summary  and  recommendations  for  further  study. 


Chapter  2 


Theory 


All  quantitative  theories  regarding  gyrotron  operation  must  consider  the  cou¬ 
pled  problems  of  the  Lorentz  force  dominated  electron  dynamics  and  the  electro¬ 
magnetic  fields  described  by  Maxwell’s  equations.  The  complete  solution  consider¬ 
ing  specific  device  geometries  is  quite  intractable,  even  using  numerical  methods. 
Therefore,  the  major  differences  in  gyrotron  theories  are  in  their  assumptions  and 
their  degrees  of  completeness.  The  linear  theory  is  useful  in  determining  the  start 
of  oscillation  conditions.  The  nonlinear  theory  is  required  to  predict  the  output 
power  of  an  oscillator,  and  therefore  is  needed  in  efficiency  calculations. 

This  chapter  will  outline  the  basic  development  of  both  the  linear  and  non¬ 
linear  theories  of  the  gyrotron.  These  theories  will  then  be  applied  specifically  to 
the  MIT  gyrotron  to  predict  its  performance  at  second  harmonic  operation.  Two 
separate  resonator  cavities  will  be  discussed.  The  first  is  described  in  Section  1.2 
of  this  paper  and  will  be  referred  to  as  the  0,3,1  cavity.  All  of  the  experimental 
work  described  in  later  chapters  used  this  resonator.  The  second  cavity  which  was 
designed  for  operation  at  the  fundamental  in  the  TE6,i,i  mode,  will  be  referred 


to  as  the  whispering  gallery  cavity.  Extensive  experimental  work  is  planned  for 
this  resonator  in  the  future.  The  chapter  concludes  by  using  a  nonlinear  theory 
to  make  an  efficiency  scaling  approximation  from  the  fundamental  to  the  second 


harmonic. 


Description  of  the  Gain  Mechanism 


The  radiation  producing  interaction  in  gvrotrons  can  be  understood  by  ele¬ 
mentary  arguments  using  a  classical  description.  The  discussion  presented  here 
follows  the  development  in  an  undergraduate  phvsics  text  by  Bekefi  and  Barrett 
33  .  The  simplified  geometry  to  be  considered  is  a  beam  of  monoenergetic  elec¬ 
trons  travelling  in  helical  paths  about  static  magnetic  field  lines.  The  electrons 
encounter  an  oscillating  electromagnetic  field  in  the  resonator  cavity,  which  for 
the  purposes  of  this  schematic  can  be  considered  unperturbed  by  the  presence  of 
the  electron  beam. 

Entering  the  interaction  region,  the  electrons  have  random  phases  in  the  plane 
perpendicular  to  the  direction  of  the  magnetic  field.  The  transverse  component  of 
the  high  frequency  RF  field  in  the  resonator  will  tend  to  accelerate  some  of  the 
electrons  and  decelerate  others.  Variations  in  the  field  across  the  Lartnor  orbit 
are  ignored  here,  so  this  discussion  is  only  applicable  at  the  fundamental.  Those 
electrons  which  are  accelerated  initially  by  the  electric  field  will  increase  in  mass 
by  the  relativistic  factor  The  increase  in  ^  will  decrease  the  cyclotron  frequency 
since  -v,  1  These  electrons  will  be  retarded  in  phase  relative  to  those  not 

initially  accelerated.  Conversely,  electrons  which  are  decelerated  by  the  electric 
field  will  have  a  resultant  decrease  in  an  increase  in  and  therefore  advance 


in  phase.  Azimuthal  electron  bunching  will  occur,  with  relativistic  mass  effects 
providing  the  bunching  mechanism. 

Simplified  conditions  for  gain  may  be  derived  i34  by  now  viewing  the  electron 
bunch  as  a  single  negative  charge  orbiting  at  u-v  near  synchronism  with  the  RF 
field,  oscillating  at  uj.  If  the  cyclotron  frequency  is  slightly  greater  than  u>,  then 
after  several  periods  of  oscillation  the  electron  bunch  will  be  positioned  such  that 
it  is  accelerating  and  gaining  energy  at  the  expense  of  the  field.  However,  if  wc  is 
slightly  less  than  then  after  several  periods  the  bunch  will  be  decelerated,  giving 
up  energy  to  the  RF  field.  In  simple  terms,  the  gyrotron  depends  on  relativistic 
electrons  for  azimuthal  bunching  and  a  gyrofrequencv  such  that  <  u.'  for  gain. 

This  description  details  the  interaction  at  the  fundamental  of  the  cyclotron 
frequency.  To  generalize  to  the  ntb  harmonic,  the  finite  Larmor  radius  must  be 
included  and  the  electron  must  be  perturbed  by  the  field  n  times  in  one  gyro  orbit. 
The  field  may  be  represented  by  an  n,h  order  multipole  expansion  and  the  beam 
interaction  involves  the  n1'1  derivative  of  the  field.  The  simplified  gain  condition 
becomes  u’  rijj, 

I  he  azimuthal  bunching  described  above  is  the  electron  cyclotron  maser  in¬ 
stability  generally  associated  with  CRM's.  Additionally,  there  is  the  possibility  of 
axial  bunching  of  the  electrons  leading  to  what  is  termed  a  Weibel-type  instability 
35  .  The  axial  bunching  arises  from  the  Lorentz  force  interaction  between  the  elec¬ 
trons  and  the  RF  magnetic  field  and  is  not  at  all  dependent  on  relativistic  effects. 
Chu  and  Hirshfield  36  have  shown  that  the  two  bunching  mechanisms  compete 
with  each  other,  so  that  depending  on  the  particular  device,  one  will  dominate.  In 


an  earlier  paper,  Hirshfield  J37]  concluded  that  the  azimuthal  bunching  dominates 
in  fast  wave  interactions.  Since  a  gyrotron  operates  near  cutoff,  the  wave  vector 
components  satisfy 


k2^  »  k?  .  (2.1) 

Here  k±  and  k\  represent  the  component  magnitudes  of  k  in  the  directions  per¬ 
pendicular  and  parallel  to  the  static  magnetic  field  respectively.  Thus  the  Doppler 
shift  term  k,,v  \  is  small  and  the  fast  wave  condition 

u.1  /  k  c  (2.2) 


is  satisfied.  Therefore  the  axial  bunching  is  ignored  in  most  gyrotron  theories. 
It  has  been  shown  ; 38 :  by  similar  arguments  that  the  electron  beam  interacts 
with  TM  modes  with  substantially  lower  efficiency  then  with  TE  modes.  For  this 
reason,  TM  modes  are  generally  not  considered  and  will  not  be  discussed  in  this 
paper. 


Definitions  of  Q  and  r)  — ■ 

It  is  useful  at  this  point  to  introduce  some  parameters  which  appear  through¬ 
out  CRM  theories  in  power  balance  relationships.  The  diffract ive  Q  (Qi>)  is  related 
to  the  stored  energy  in  the  cavity  (/')  and  the  RF  output  power  (/J,„lt)  by 


Qn 


(2.3) 


Since  some  energy  is  extracted  by  ohmic  losses,  a  similar  relationship  exists  be¬ 
tween  the  ohmic  Q  (Qohm)  and  the  power  that  is  used  in  wall  heating  (F,.iin,)- 
The  overall  quality  factor  (Q)  includes  both  diffractive  and  ohmic  losses. 


Q  Qn  '  Q  .u» 
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The  output  efficiency  77  of  a  gyrotron  is  defined  by  77  =  P0ut/P\n  where  P\n 
represents  the  beam  power  given  by  Pin  —  IV c.  Here  /  is  the  beam  current  and  V'c 
the  cathode  voltage.  As  with  the  Q  factors,  consideration  of  the  fact  that  some 
of  the  RF  energy  heats  the  walls  instead  of  travelling  out  the  waveguide  leads 
to  the  definition  of  the  electronic  efficiency  rjei >.  Both  power  loss  mechanisms  are 
included  in  r]ek  by  the  expression  Tjet  —  (/Cut  1  Pobm) / P\n-  The  output  efficiency 
may  be  written  in  terms  of  rjet  by 


Since  only  TE  modes  will  be  considered  and  the  static  magnetic  field  can  do  no 
work  on  the  electrons,  only  the  perpendicular  velocity  (uj  )  of  the  electrons  can 
contribute  energy  to  the  field.  The  parallel  electron  velocity  (r,  )  remains  constant 
throughout  the  cavity  to  the  extent  that  the  RF  magnetic  field  may  be  ignored. 
The  final  efficiency  to  be  defined,  77^,  accounts  for  the  fad  that  only  the  transverse 
beam  energy  may  be  converted  to  electromagnetic  field  energy.  Equation  (2.0) 
expresses  the  relationship  between  7/,,  and  7/  . 


1 


(2.0) 


As  a  final  set  of  useful  definitions,  the  normalized  velocities  are  3  v _  c 
and  3  v  c.  They  are  related  to  3  introduced  in  Chapter  1  by  3 2  ~  3~. 

The  relativistic  y  written  in  terms  of  3  in  equation  (1.3)  may  be  calculated  by 


1  -  1  - 


V,. 

511 


(2.7) 


when  \\  is  in  kilovolts. 


2.1  —  The  Electromagnetic  Field  Structure  in  an  Open  Resonator 


All  of  the  theories  to  be  considered  require  some  knowledge  of  the  RF  field 
structure  in  the  cavity.  The  complete  solution  of  this  problem  in  actual  gvrotron 
operation  involves  a  coupled  system  of  equations  relating  the  perturbing  effect 
of  the  electron  beam  on  the  cavity  fields  and  the  RF  field  effects  on  the  particle 
dynamics.  A  reasonable  first  approximation  is  the  cold  cavity  electromagnetic  field 
solution  disregarding  the  presence  of  the  electron  beam.  The  CAVRF  computer 
code  39  developed  at  the  Naval  Research  Laboratory  provides  such  a  solution. 

The  gyrotron  resonator  is  basically  a  section  of  irregular  waveguide  aligned 
with  the  magnetic  field.  The  solution  to  Maxwell’s  equations  without  source  terms 
in  a  perfectly  conducting,  nonuniform  waveguide  had  been  considered  previously, 
but  Vlasov  et  al.  -10  provided  an  approach  specifically  applicable  to  gvrotrons. 
In  39  .  Fliflet  and  Read  continue  with  this  development,  rederiving  the  applicable 
equations  and  discussing  the  numerical  methods  employed  in  creating  the  CAVRF' 
program. 

In  outlining  the  CAVRF  solution  technique,  the  longitudinal  axis  parallel  to 
the  magnetic  field  will  be  ifie  '-direction.  The  longitudinal  profile  function  of  the 
field  /...(.;)  satisfies  the  one  dimensional  wave  equation 

(ft  f 

(,'2  -  k;(^.=)L  0  .  (2.8) 

In  equation  (2.8).  k..  is  the  complex  equivalent  to  k  and  the1  “s"  subscripts  denote 
t  hat  the  quantities  are  mode  dependent .  The  radiation  boundary  conditions  which 


ensure  that  the  fields  are  decaying  exponentially  in  the  cutoff  input  section  of  the 
cavity  and  are  travelling  waves  in  the  output  section  are 


and 


(2.9) 


Here  2jn  and  2oll,  refer  to  the  coordinate  positions  of  the  input  and  output  ends 
of  the  cavity. 

The  solutions  to  equation  (2.8)  with  the  boundary  conditions  (2.9)  yield  com¬ 
plex  eigenfrequencies.  The  real  part  of  the  eigenfrequencv  is  related  to  the  actual 
frequency  of  radiation  ( i /)  in  the  standard  manner. 


Re  u.* 
v  — - 
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(2.10) 


The  imaginary  part  can  be  used  to  determine  the  diffractive  Q  of  the  resonator 
by  the  relation 


Qd  = 


Re 


(2.11) 


2(Im  w) 

The  CAYRF  program  separates  equation  (2.8)  into  real  and  imaginary  parts 
and  numericailx  integrates  the  resulting  pair  of  coupled  second  order  differential 
equal  ion-;.  The  eigenfrequencv  is  found  by  an  iterative  procedure  based  on  the 
method  of  steepest  descent.  The  code  requires  the  actual  cavity  dimensions  as 
inputs  and  yields  frequency,  diffractive  Q.  and  the  magnitude  and  phase  of  f?  as 
a  function  of  .:. 


Application  of  the  CAVRF  Code 


Figure  2.1a  shows  a  plot  of  the  normalized  field  magnitude  as  a  function 
of  cavity  length  using  an  actual  CAYRF  output.  The  result  depicted  is  for  the 


cavity  length  (CM. 


TE 1 1,2,1  mode  using  the  0,3.1  cavity  dimensions.  Since  this  cavity  was  designed 
to  oscillate  at  140  GHz  in  the  TEu,3.i  mode,  it  has  a  frequency  of  241  GHz  in  the 
TEn^.i  mode.  Superimposed  over  the  field  magnitude  plot  are  straight  line  seg¬ 
ments  representing  the  0,3,1  cavity  radius  as  a  function  of  z.  The  radial  dimension 
scale  on  the  right  is  not  the  same  as  the  longitudinal  scale  on  the  horizontal  axis. 

The  CAVRF  code  remains  sufficiently  accurate  only  for  resonator  cavities 
which  have  a  smoothly  varying  radius  along  c.  Abrupt  changes  in  the  cavity  radius 
violate  the  "weakly  irregular"  condition  which  permits  separation  of  variables  in 
solving  the  Helmholtz  equation.  Both  the  0.3.1  and  whispering  gallery  cavity 
profiles  should  satisfy  the  smoothly  varying  stipulation.  The  downtaper  on  the 
input  end  is  0.15  for  the  0.3.1  cavity  and  0.5  for  the  whispering  gallery  cavity. 
Both  resonators  have  output  uptaper  angles  of  4  . 

The  data  presented  in  Tables  2.1a  and  2.1b  was  generated  by  the  CAVRF 
code.  Table  2.1a  represents  all  the  possible  TEm,;,,i  modes  which  could  oscillate 
at  the  second  harmonic  between  4.0  and  G.O  tesla  in  the  0.3.1  cavit\.  Table  2.11) 
includes  all  the  second  harmonic  TE„i  ;  .i  modes  in  the  whispering  gallery  cavity 
o\er  the  same  magnetic  field  range.  Using  the  TEm-;.w  notation,  only  q  1  modes 
were  considered.  Temkin  41  has  shown  that  Qi,  ~~  q  2  for  high  Qp  cavities,  so 
higher  axial  modes  (q  1 )  would  have  substantially  reduced  diffractive  Q's.  It  will 
be  evident  later  in  this  chapter  that  these  lower  Qp  modes  will  require  significantly 
higher  beam  currents  for  the  start  of  oscillation  condition  and  w’ill  produce  much 
lower  efficiencies.  In  the  overmoded  cavities  under  consideration,  these  higher 
order  q  modes  should  nol  be  present. 


Table  2.1a 

—  0,3,1  Cavity  u  ss  2u >c 

CAVRF  Results 

MODE 

FREQUENCY  (GHz) 

Qd 

*!!  (cm  ') 

TE,  5  i 

203.72 

3310 

1.537 

TEi3,l,] 

204.61 

3350 

1.538 

TE63.i 

209.26 

3510 

1.545 

TE9,2,i 

209.51 

3520 

1.545 

TE4  4,1 

218.78 

3870 

1.559 

TE,41  1 

218.94 

3880 

1.559 

TEj  5,i 

224.03 

4080 

1.566 

TE10,2.i 

225.40 

4140 

1.568 

TEo.5.1 

225.71 

4150 

1  568 

TE7.3] 

226.52 

4180 

1.569 

TE  is,]  ,i 

233.24 

4460 

1.578 

TEb.4,1 

237.24 

4630 

1.583 

TE]  12  1 

241.17 

4800 

1.587 

TE8,3.1 

243.55 

4910 

1.590 

TE35 1 

243.75 

4920 

1.591 

TEi  6.1 

246.86 

5060 

1 .594 

TEi6.ii 

247.51 

5090 

1.595 

TE6  4,1 

255.37 

5450 

1.603 

TE12  2  1 

256.84 

5530 

1.605 

TE9.3  1 

260.39 

5700 

1 .609 

TEit.i.i 

261.76 

5770 

1.610 

TE4.fi.! 

263.01 

5830 

1.611 

TE2  6  1 

267.35 

6040 

1.616 

TE„.6  1 

268.76 

6120 

1.617 

TEjs.2.1 

272.12 

6310 

1.621 

TE7  1  1 

273.22 

6360 

1 .622 

TEi*.,  1 

275.99 

6490 

1 .624 

TEjo  3  1 

277.07 

6660 

1 .625 

TEr,  fi  1 

281.89 

6820 

1 .630 

TE3.6.J 

287.33 

7120 

1 .635 

TE14.2 . 

287.92 

7140 

1 .635 

TE,  7., 

289.95 

7260 

1 .637 

TE19  1  , 

290.20 

7280 

1 .637 

TE8  4  i 

290.84 

7320 

1 .638 

TE, 1,3.1 

293.60 

7470 

1 .640 

TEe  5.1 

300.46 

7680 

1 .652 

TEib  2,1 

303.35 

7690 

1.660 

30 
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The  possible  modes  for  the  given  magnetic  field  range  were  found  by  scaling 
ump  from  a  fundamental  mode  where  i/mp  represents  the  pth  root  of  J'm{x)  =  0. 
Here  J'm{x )  is  the  derivative  with  respect  to  argument  of  the  Bessel  function  of 
order  m.  The  cavity  dimensions  required  as  CAVRF  input  values  could  not  be 
measured  to  the  desired  accuracy,  so  the  CAVRF  output  frequencies  for  the  second 
harmonic  modes  were  scaled  by  the  ratio  of  experimentally  measured  frequency 
to  CAVRF  frequency  of  a  known  fundamental  mode.  The  accuracy  of  the  ex¬ 
perimental  technique  discussed  in  Chapter  3  allowed  five  significant  digits  in  the 
scaled  CAVRF  frequencies  tabulated  in  Tables  2.1a  and  2.1b.  Since  Qp  was  not 
determined  experimentally  at  the  fundamental  with  sufficient  accuracy,  the  direct 
CAVRF  output  values  of  Qp  at  the  second  harmonic  are  presented  in  the  tables. 

--  Using  CAVRF  for  the  Gaussian  Approximation  to  f(z)  — 

The  final  columns  in  Tables  2.1a  and  2.1b  are  labelled  k\ which  represents  a 
parameter  required  for  theories  discussed  laler  in  this  chapter.  The  field  magni¬ 
tude  profile  plotted  in  Figure  2.1  resembles  a  Gaussian  function  with  some  degree 
of  imagination.  This  Gaussian  approximation  is  used  in  several  theories  In  sub- 
si  it  lit  ing  t  lie  function 

f(=)  e  l2:  L  "y  (2.12) 

for  the  actual  cold  cavity  /..(c)  solved  for  In  CAVRF.  Here-  lJV g  represent s  the 
effective  cavity  length,  or  the  full  width  of  the  Gaussian  function  at  1  he  1  c  point. 
This  is  related  to  the  tabulated  value  of  k  b\ 


The  mechanics  of  calculating  the  k  values  involved  running  CAVRF  with 
975  data  points  along  the  2-axis  for  sufficient  accuracy.  The  highest  three  field 
magnitude  values  and  corresponding  2-coordinates  were  locally  fit  to  a  parabola  to 
determine  the  maximum  point  and  its  2-coordinate  (2max)-  A  parabola  was  then 
fit  to  the  three  points  closest  to  the  1  /e  value  of  the  maximum  on  the  input  side, 
and  the  2-coordinate  (2]  ,e)  was  determined.  The  input  side  was  used  because  the 
efficiency  is  known  to  depend  strongly  on  the  RF  field  profile  on  t  he  input  side  and 
only  weakly  on  the  profile  on  the  output  side.  The  k  value  was  then  calculated 
by 

k  -  -—  1  --  .  (2.14) 

This  value  is  close  but  not  necessarily  the  same  as  the  parallel  wave  number 
k  expressed  as 

k  -  y  (*•  t)2  -  k'i  .  (2.15) 

In  typical  cylindrical  resonator  geometries. 

k.  u,„v  r,.  (2.10) 

w  here  r  .  is  the  cavity  radius.  This  implies  k_  is  a  function  of  2.  As  shown  in 
Figure  2.1.  the  fields  are  generally  confined  to  the  straight,  central  region  of  the 
cavity.  In  practical  calculations  of  k  .  it  suffices  to  use  t lie  straight  section  radius 
in  equation  (2.10).  The  condition  k~  k2  implies  (.*.-  c)*  -  k2  k 2  .  For 


this  reason,  equation  (2.15)  is  not  very  useful  uidess  (.2  r)  and  k  are  known  to 
unrealistic  accuracies.  The  parameter  k  described  above  actually  provides  a  more 


reasonable  estimation  of  k\i.  Unless  stated  otherwise  in  this  paper,  hi,  will  be  used 
whenever  a  value  of  k t  is  required. 

The  CAVRF  results  provide  an  adequate  starting  point  estimation  of  the 
electromagnetic  fields  in  the  resonator.  However,  gyrotron  performance  theories 
require  consideration  of  the  energy  exchange  between  the  electrons  and  the  fields. 
The  linear  theory  will  serve  as  a  first  step  in  this  regard. 

2.2  —  The  Linear  Theory  of  the  Gyrotron 

To  define  the  linear  regime  of  CRM  operation,  the  electron  beam  current 
density  ./  is  in  general  not  only  the  unperturbed  beam  density  Jt..  but  a  sum  of 
terms  describing  the  coupling  between  the  beam  and  electromagnetic  fields.  Thai 
is.  ./  .7  -  C]£(<f)  *  coS~(S)  -  c ?,f3(f)  -  •••  where  f  is  the  field  amplitude 

and  f  a  unit  vector  in  the  direction  of  f .  For  small  field  values,  the  second  order 
and  higher  terms  may  be  neglected,  reducing  J  to  terms  linear  in  f .  Using  the 
output  power  relation  P,,u t  ~  J  J  ■  f ,  it  follows  that  in  the  linear  limit,  t;  if*'. 

The  small  field  amplitude  condition  restricts  the  range  of  applicability  to  the 
threshold  of  mode  oscillation.  The  primary  application  of  the  linear  theory  is 
the  <  al<  illation  of  the  starting  current  /  T .  which  is  the  minimum  beam  current 
required  for  mode  initiation.  This  allows  the  determination  of  which  modes  should 
be  excited  at  a  specific  set  of  experimental  operating  conditions  including  beam 
current  /.  cathode  voltage  V,.,  and  magnetic  field  value  B.  The  utility  of  the  linear 
theory  is  enhanced  because  it  generally  results  in  analytic  expressions,  as  opposed 
1o  the  nonlinear  theory  which  must  bo  solved  numerically. 


There  have  been  a  number  of  formulations  of  the  linear  theory,  differing 
mainly  in  their  scope  of  applicability.  In  1965,  Hirshfield  |37  developed  a  the¬ 
ory  for  the  TEn.1,1  mode  at  the  fundamental.  Chu’s  treatment  42!,  restricted  to 
the  azimuthally  symmetric  TEo1?,.<j  modes,  extends  the  theory  to  cyclotron  har¬ 
monics.  Temkin  et  al.  ;43j  derived  an  expression  for  /ST  at  the  fundamental  for 
arbitrary  TEm,?.w  modes.  All  of  these  methods  assumed  a  sinusoidal  longitudinal 
distribution  of  the  RF  field.  The  most  comprehensive  Soviet  investigation  was 
performed  by  Petelin  and  Yulpatov  44 1,  but  a  particularly  simple  yet  meaningful 
expression  was  reported  by  Gaponov  45  which  will  be  used  extensively  later  in 
this  section.  The  Gaponov  equation  assumes  a  Gaussian  field  distribution  and 
yields  / . T  at  the  magnetic  field  value  corresponding  to  the  minimum  1<T. 

This  section  will  outline  the  development  of  the  linear  theory  as  formulated 
by  Kreischer  and  Temkin  46!.  Their  treatment  encompasses  all  TE  modes  in¬ 
cluding  harmonics,  and  any  longitudinal  field  profile.  Additionally,  an  arbitrary 
distribution  of  electron  velocities  may  be  included.  In  the  appropriate  limits, 
their  expressions  reduce  to  those  results  cited  earlier.  This  will  be  demonstrated 
explic  itly  in  t  he  case  of  the  Gaponov  equation. 

Description  of  the  Line  r  Analysis 

Several  assumptions  are  necessary  in  the  Kreischer  formulation  of  the  linear 
theory.  These  include  the  same  condition  required  by  CAVRF  that  the  cavity 
radius  be  slowly  varying  along  the  r-axis.  The  electron  beam  should  be  weakly 
relativistic,  since  the  theory  only  retains  the  relativistic  factor  q  in  the  denominator 
of  Space  charge  effects  are  neglected  as  is  the  interaction  between  the  electrons 


and  the  RF  magnetic  field.  This  latter  statement  eliminates  consideration  of  axial 
bunching  and  TM  modes  as  discussed  previously. 


The  derivation  is  based  on  the  combined  Vlasov  and  Maxwell  equations 
as  a  foundation  for  describing  the  beam  and  electromagnetic  field  interaction. 
Maxwell's  equations  are  solved  by  a  Slater  expansion  of  the  fields  in  terms  of  or¬ 
thonormal  modes.  When  the  cavity  is  at  equilibrium,  the  power  balance  equation 


mav  be  written 


u,’s  11 

-  -  Re  n 

Q .  fl,  P. 


(2-17) 


where  P,  represents  the  time  independent  Sialer  expansion  coefficient  for  the  mode 
The  “.s"  subscripts  on  Q.  and  £  indicate  that  they  are  also  mode  dependent. 
This  equation  follows  from  the  total  Q  generalization  of  equation  (2.3).  In  this 
case  the  stored  energy  l'  =  / 2)  and  the  total  power  extracted  from  the  beam 

is  the  interaction  integral  f  J  ■  £f.  Since  (2.17)  is  applicable  at  equilibrium,  il 
is  representative  of  the  oscillation  threshold  conditions.  Solving  the  equation  for 
the  current  yields  1  T. 

The  Vlasov  equation  is  linearized  in  terms  of  the  fields  and  electron  distribu¬ 
tion  function.  Tint  zero  order  equation  b  subtracted  from  the  perturbed  equation, 
and  the  resulting  expression  is  solved  for  the  perlurbed  electron  distribution  func¬ 
tion  / 1  by  the  method  of  characteristics.  The  J  in  equation  (2.17)  may  then  be 


solved  for  by  the  equation 


J  =-  e 


I  dv  vN  f]  (rj  . 


(2.18) 


The  integration  in  equation  (2.18)  is  over  velocity  and  N  is  the  spatial  beam 
density.  The  general  interaction  integral  that  results  from  these  calculations  is 
quite  complicated,  involving  several  geometric  factors  and  functions  which  are 
determined  by  the  longitudinal  field  profile.  The  general  expression  will  not  be 
presented  here. 

The  derivation  at  the  fundamental  assumes  thal  the  electron  interacts  with 
the  field  value  at  the  gyrocenter.  To  extend  the  theory  to  include  higher  harmon¬ 
ics.  the  field  value  on  the  Larmor  orbit  must  be  used.  To  accomplish  this,  the 
scalar  potential  is  expanded  in  terms  of  four  coordinates.  These  are  the  cylin¬ 
drical  coordinates  azimuthal  angle  and  radius  (rf)  of  the  guiding  center,  as  well 
as  the  phase  and  Larmor  radius  (r,.)  of  the  electron  about  the  gyrocenter.  The 
transverse  field  is  then  determined  from  the  scalar  potential  and  the  resi  of  the 
derivation  remains  the  same  as  for  the  fundamental.  Again,  the  general  case  will 
not  be  shown,  but  instead  the  specialized  result  for  an  infinitely  thin,  azimuthally 
symmetric,  annular  electron  beam  with  radius  r,  is  presented. 


(..  ~  p  m'w.  k~v  (dJ„{  k  .  rL)  \ 
2~  (J  c  Cp  l  d(k,rL)  ) 


/  -  , 

1  /  rJ  \  <IF, 
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1  dr  .  /,,  | 

(  r  ) 

F 

2  \  A'  v  c  2  /  dx 

In  equation  (2.19).  /.  is  the  unperturbed  electron  distribution,  r  is  the  average 
parallel  velocity  and  all  of  the  mode  subscripts  have  been  dropped.  Thefunclion  F, 
depends  on  the  longitudinal  field  profile  function  f(z)  and  may  be  calculated  in 
terms  of  the  Fourier  transform  of  /  with  the  independent  variable  .r.  This  .r  is 


known  as  the  detuning  parameter  and  is  a  measure  of  the  off-resonance  condition. 


Explicitly,  x  is  defined  by 


TtsjJ  r' 


(2.20) 


and  it  may  be  considered  a  normalized  value  of  B. 

The  G d  factor  measures  the  effectiveness  of  the  interaction  between  the  elec¬ 
tron  beam  and  the  fields.  It  is  evaluated  as  G p  J;n^Tt(k±re).  The  “±”  sign 
is  necessary  because  there  is  a  double  degeneracy  in  the-  modes  with  azimuthal 
mode  index  m  0.  This  corresponds  1o  two  senses  of  azimuthal  rotation  of  the 
mode.  There  is  a  lack  of  symmetry  in  the  problem  due  to  the  specific  direction 
of  rotation  of  the  electrons  about  the  magnetic  field  lines.  Correspondingly,  the 
Hesse]  function  values  of  J?n^n(k^r£)  and  J%.t_ri(k_  rf)  may  differ  greally  yield¬ 
ing  significantly  different  values  of  Lr.  The  degeneracy  is  not  resolved  in  the 
azimuthally  symmetric  (rrt  —  0)  modes  which  are  standing  waves  in  the  azimuthal 
direction.  Mathematically  this  is  realized  by  J?n{k^r()  ~  J‘„(k_  rt).  thus  yield¬ 
ing  the  same  /.T  regardless  of  sign. 

Derivation  of  the  Gajtonov  Result 

Am  arbitrary  distribution  of  velocities  max  !»<■  treated  in  equation  (2.1!))  In 
specifying  /,  .  In  the  case  of  no  velocity  spread,  f,  becomes  the  delta  function 
/•  *(»\  r  „)<“)()'  r  .,).  and  r  ?•  ,  .  The  integral  in  the  denominator 

of  equation  (2.19)  may  now  be  evaluated. 
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(2.21) 


as 


Equation  (2.21)  is  now  appropriate  for  a  gyrotron  magnetron  injecton  gun  produc¬ 
ing  a  perfectly  thin,  annular  electron  beam  with  no  velocity  spread.  To  proceed 
further,  the  longitudinal  field  distribution  must  be  specified.  By  assuming  the 
Gaussian  function  for  f(z)  from  equation  (2.12)  and  evaluating  the  expression 
above  at  the  x  value  corresponding  to  the  minimum  /.T,  the  Gaponov  equation 
45  will  result. 

As  shown  in  46  ,  the  stored  energy  P  2  for  a  Gaussian  profile  is 

?2=  (\Y  'y'U  -  J?nWmV)  •  (2.22) 

Specifxing  a  Gaussian  also  determines  F,  and  its  derivative  with  respect  to  x. 
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(2.23) 


Dropping  the  V’  subscripts  on  r  ,  and  ,  and  substituting  in  expressions  (2.22) 
and  (2.23).  the  equation  for  V.T  becomes 
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/  L^k2r2  \  u;ru  m-  Jj(iynq  ) 
\  "**:  "  /  n„n(k\rf) 


(2.24) 


I  in  approximations  applicable  to  gyrotrons.  k2  sr  k2  and  u.'  fts  nix.',.,  permit  the 
substitution  k  r,  n.i  .  The  small  argument  approximation  for  the  Bessel 
function  may  lie  used  in  the  weakly  relativistic  limit  to  yield 

/  d.l„(k  rL)  ^  /  d.Jn(r>;1  )  mull  I  /  n"  \  ‘  /  ,4  \  2(''  1 
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As  discussed  in  Section  2.1,  for  the  Gaussian  k\<  =  At||  =  Lefr/2.  This  relationship 
with  the  near  cutoff  and  near  resonance  approximations  allow 


Prior  to  making  these  substitutions,  the  normalized  cavity  length  //  and  the  beam 
interaction  term  G  will  be  introduced.  These  definitions  appear  in  numerous  Soviet 
papers  including  <8.  and  45  . 


'?2  (  Left 


H  -Z  IT 


3 


(2.27; 


G  =  - 


Jm±rAk,rt) 


vmV  ~  m2' 


(2.28) 


Using  equations  (2.25)  -  (2.28),  the  expression  for  I< T  becomes 
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(2.29) 


x// 


The  detuning  value  which  produces  the  minimum  starting  current  (xM,j,;)  is 
determined  In  taking  the  derivative  of  equation  (2.29)  with  respect  to  x  and  setting 
it  equal  to  zero.  The  resulting  value  is  given  by 


(*)  +  V  (")' 
V  W  \  \  e  / 


1 


(2.30) 


Only  minor  cosmetic  changes  remain  to  convert  (2.29)  to  the  exact  form  of 
the  Gaponov  equation.  These  include  defining  “x  Gaponov"  (x0)  as 


Tc  — 
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v  ^  /  V  v M  > 


(2.31) 


and  multiplying  numerator  and  denominator  by  10~3.  Evaluating  the  numerical 
factors  and  physical  constants  in  MRS  units  yields  2.16  amps  when  multiplied  by 
the  10“ 3.  It  appears  that  Gaponov  combined  the  2.16  and  -y  to  get  2.3,  which 
is  accurate  at  Vc  ss  32  kV.  The  following  equation  is  the  Gaponov  result  with 
definitions  (2.27),  (2.28)  and  (2.31). 
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Application  of  the  Linear  Theory  — 

The  denominator  of  the  expression  above  contains  the  total  Q ,  which  may  be 
calculated  by  equation  (2.4).  The  diffractive  Q  is  obtained  from  the  CAVRF  results 
of  Tables  2.1a  and  2.1b.  Equation  (2.33)  below  yields  Q0hm  for  the  cylindrical 
geometry  resonators  under  consideration. 
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(2.33) 


Here  r,  is  the  radius  and  />  the  skin  depth  of  the  conducting  cavity  wall.  To 
present  order  of  magnitude  examples  of  Q>. .j, ,,,  for  the  0.3.1  cavity,  the  "volume" 
modes  (radial  mode  index  p  >  3)  all  had  Q,.hm  >  20.000.  The  TE„,.^.|  modes 
with  w  >  p  -  2  had  Q,,hni  =  16,000  and  the  whispering  gallery  modes  (TE,„.t.i 
with  m  .y  1)  had  Q.,inn  =:  6.000.  These  values  of  Qohm  combined  with  the  values 
of  Qi,  ranging  from  3310  7680,  imply  that  with  the  possible  exception  of  the 
whispering  gallery  modes  with  p  -  1,  Q  ss  Qp.  This  serves  to  justify  the  previous 
neglect  of  higher  order  axial  modes  since  /  T  ~  Q  1  ~  q2.  With  the  given 


range  of  Q  values,  TEm)Pi2  should  have  approximately  four  times  the  starting 
current  as  TEmi?7,i. 

All  of  the  factors  appearing  in  equation  (2.32)  may  now  be  calculated  by 
the  data  produced  from  CAVRF  and  the  specific  experimental  conditions.  The 
frequency  and  fcy  from  Tables  2.1a  and  2.1b  yield  Leff  and  A.  An  operating  voltage 
of  Vc  —  64.3  kV  and  the  MIT  electron  gun  design  value  of  {Q1j0[)  =  1.49, 
produces  3±  =  0.383  and  /?,|  =  0.256.  These  values  permit  the  calculat  ion  of  p  by 
equation  (2.27)  and  x-  by  specifying  the  harmonic  number  n  2  in  equation  (2.31 ). 

The  frequency  from  CAVRF  also  appears  indirectly  in  the  calculation  of  Q,  by 

1 

the  v  -  dependence  of  the  skin  depth  appearing  in  equation  (2.33)  for  Q.  hm- 

The  remaining  factor  in  the  Gaponov  equation  (2.32)  is  the  beam  interaction 
parameter  G  defined  by  (2.28).  The  Ay  in  the  argument  of  the  Bessel  function 
in  the  numerator  of  (2.28)  is  determined  by  equation  (2.16).  The  specific  wall 
radius  values  used  were  r,,  -  0.348  cm  for  the  0.3.1  cavity  and  r0  =  0.250  cm 
for  the  whispering  gallery  cavity.  Both  cavity  calculations  used  a  beam  guiding 
center  radius  of  r,  -  0.197  cm.  The  degeneracy  in  G  discussed  earlier  for  the  non- 
azimuthallv  symmetric  modes  requires  an  additional  term  in  ihe  inode  notation. 
In  this  paper,  TEm<?N1?(+)  refers  to  the  positive  sense  of  rotation  and  therefore  a  G 
calculated  with  J%l^n(k±re),  while  TEmw,  ^(._ ^  is  the  opposite  with  G  resulting 


from  J'f):  „(A-  rj.  The  azimuthally  symmetric  modes  will  remain  TE,1<?,W. 

Starting  currents  were  calculated  by  equation  (2.32)  for  all  the  modes  listed 
in  Tables  2.1a  and  2.1b.  The  results  are  depicted  in  graphical  form  in  Figures  2.2a 
and  2.2b  and  the  actual  values  are  lisled  in  Tables  2.2a  and  2.2b. 
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Second  Harmonic  Starting  Currents  for  the  0,3,1  Cavity 


MAGNETIC  FIELD  (tesla) 


Table  2.2a  —  0,3,1  Cavity  a >  «  2u >c  Starting  Currents 


MODE 

B  (tesla) 

Q 

Eeff  /  A 

1st  fam; 

te1<BiJ(+) 

4.054 

2910 

8.845 

1.4 

TE13, ,.!(-) 

4.072 

2120 

8.876 

41. 

TE6i3.,(+) 

4.165 

3000 

9.035 

1.1 

TEfl^.jf-) 

4.170 

2880 

9.044 

0.72 

TE4i4ii(4) 

4.356 

3320 

9.364 

2.1 

TEl4,l,l(_) 

4.360 

2312 

9.369 

83. 

TEa,B,,(  +  ) 

4.462 

3500 

9.545 

1.1 

TEjo,2.1(-) 

4.489 

3280 

9.592 

0.66 

TEo  5.1 

4.496 

3560 

9.603 

1.9 

TE7.3.H-r) 

4.512 

3480 

9.630 

1.2 

TEjs.j  k  ) 

4.647 

2500 

9.861 

110. 

TE5  4.1(^  1 

4.727 

3870 

10.00 

1.1 

TEj j  2.i(  ) 

4.806 

3680 

10.14 

0.66 

TE8  3  ,it_  j 

4.854 

3970 

10.22 

1.3 

TE3lB.|(4) 

4.858 

4120 

10.22 

1.1 

TEa>6>1(+) 

4.921 

4240 

10.33 

5.6 

TEa6, !,!(_) 

4.934 

2680 

10.35 

180. 

te6.4i1(+) 

5.092 

4440 

10.62 

0.68 

TE12.2,1(-  ) 

5.121 

4090 

JO. 68 

0.71 

TE9  3  !(_,.) 

5.193 

4480 

10.80 

1 .5 

TE17,  l.l(-) 

5.220 

2850 

10.85 

300. 

TE4.5.J(j  j 

5.245 

4760 

10.89 

1.5 

TE2,o,](  ) 

5.332 

4910 

11.04 

1.8 

TEo.r,  1 

5.361 

5000 

1 1 .09 

3.5 

TE  13.2.  Ji  ) 

5.434 

4  500 

11.21 

0,81 

TE7  4  l(, ) 

5.450 

5030 

1 1 .24 

<),5i  1 

TEjs.i.k _ ) 

5.506 

3000 

11.34 

660 

TE10.3.k  .  ) 

5.528 

5000 

11.37 

1.7 

TE5.5 .!(_  I 

5.625 

5420 

11.54 

0.78 

TE3.6.K  -  ) 

5.734 

5670 

11.73 

0,75 

TEj4  2.1(  ) 

5.746 

4910 

11.75 

0.97 

te1iT.1(  , 

5.787 

5780 

1 1 .82 

0.77 

TEio.j.ii  ) 

5.792 

3150 

11.83 

890. 

TEs.4.j(^  ) 

5.804 

5630 

11.85 

0.41 

TEj  1 .3  1  (  ) 

5.860 

5510 

11.91 

0.76 

TEfi  5,]| _) 

5.998 

5970 

12.13 

0.72 

Table  2.2b  —  Whispering  Gallery  Cavity  u;  ss  2u >c  Starting  Currents 


MODE 

B  (tesla) 

Q 

Leff/\ 

1st  (am 

TE9.1K- ) 

4.029 

606 

5.887 

4.0 

TE3  3 k  _ ) 

4.270 

722 

6.118 

11. 

TEj  4  1(_  | 

4.407 

764 

6.244 

15. 

TE6  2.1(  4  ) 

4.418 

751 

6.257 

6.1 

TEio.ui- ) 

4.431 

693 

6.266 

3.6 

TE4  3.1  ( - ) 

4.777 

863 

6.596 

8.6 

TEn  j ,i(_) 

4.832 

780 

6.648 

3.4 

TE7i2.1(  +  ) 

4.873 

875 

6.690 

5.3 

te2.4.1(-) 

4.963 

915 

6.776 

9.0 

TEo.4,1 

5.021 

933 

6.832 

11. 

TE12,i,i(-) 

5.232 

866 

7.032 

3.3 

TE6,3,l(_) 

5.274 

1000 

7.075 

8.4 

TE8,2i1(  +  ) 

5.322 

1000 

7.120 

4.7 

TE3i4.1(-) 

5.501 

1080 

7.296 

8.1 

TEj,5,1(  ) 

5.607 

1110 

7.396 

™  r 

i  .«> 

TE]  3  j  1  (  ) 

5.632 

952 

7.420 

•)  n 

TE0.3.1 1  -  ) 

5.761 

1150 

7.550 

10 

TEgj.K-.  ) 

5.768 

1120 

7.556 

4.3 

In  the  figures,  each  vertical  line  represents  an  individual  mode  plotted  at  the  B 
field  value  corresponding  to  xm;n.  Equation  (2.30)  determined  xmjn,  and  for  all 
modes  the  value  of  xn,in  ~  —1.1.  This  is  a  result  of  the  one  under  the  square 
root  sign  being  the  dominant  term  since  n  »  n.  The  xmin  value  is  converted 
to  B  by  equations  (2.20)  and  (1.2).  The  vertical  axis  has  an  inverted  log  scale 
so  the  modes  with  the  lowest  starting  currents  appear  as  the  longest  lines.  The 
horizontal  dashed  line  represents  the  experimental  operating  current  of  5  amps. 
Modes  extending  above  this  line  should  be  above  threshold.  Space  did  not  permit 
labelling  each  mode,  but  correlating  the  B  and  7.T  values  in  the  accompanying 
tables  allows  individual  mode  identification. 

The  values  presented  in  Tables  2.2a  and  2.2b  include  the  two  parameters  Q 
and  Letf  A  which  were  used  in  equation  (2.32).  The  accuracy  in  the  starting 
current  calculation  limits  the  tabulated  values  of  JST  to  two  signigicant.  digits.  Only 
one  sense  of  mode  rotation  is  listed  in  the  tables.  Values  of  />T  were  computed 
for  both  “(  — )”  and  “(-)*  modes,  the  result  appearing  in  the  tables  representing 
the  mode  with  the  lowest  LT.  Generally  the  two  calculations  yielded  values  which 
differed  by  at  least  an  order  of  magnitude.  The  mode  with  the  higher  1  T  would  be 
in  the  excitation  range  of  the  mode  with  the  lower  starting  current,  and  therefore 
the  higher  /.T  mode  should  not  be  present.  The  0,3,1  cavity  TE],?,,]  was  the  only 
mode  with  the  same  “(  +  )”  and  “(-)”  I<T  values  to  two  significant  figures.  The 
fourth  digit  resolved  the  issue  in  favor  of  the  TEj  +  )  mode,  but  it  would  not 
be  considered  a  reliable  distinction. 


Examination  of  Figure  2.2a  reveals  the  overmoded  nature  of  the  0,3,1  cavity. 
There  are  36  modes  with  B(xmin)  values  between  4.0  and  6.0  tesla.  Of  these, 
28  modes  are  predicted  by  the  Gaponov  equation  to  be  above  threshold.  The 
exceptions  are  the  whispering  gallery  modes,  with  starting  currents  ranging  from 
41  to  890  amps,  and  the  TE1j6j1(+)  with  /ST  =  5.6  amps.  There  are  14  modes 
listed  with  a  predicted  starting  current  of  less  than  an  amp.  This  would  point  to 
the  possibility  of  TEm,fl.2  modes  by  the  simple  scaling  of  Q  values,  but  the  mode 
density  is  prohibitive. 

By  contrast,  the  whispering  gallery  cavity  spectrum  is  fairly  sparse  as  seen  in 
Figure  2.2b.  This  was  a  design  consideration  to  reduce  mode  competition  at  the 
fundamental.  The  reduced  wall  radius  limits  the  number  of  possible  modes  to  18 
in  the  4.0  and  6.0  tesla  range.  Only  seven  modes  appear  above  threshold  on  the 
plot,  those  being  the  whispering  gallery  modes  and  two  other  “surface”  modes, 
the  TE8.2.i(  +  )  and  the  TEp,2.i(+).  The  range  of  ]<T  values  is  much  more  limited 
in  this  case,  with  the  lowest  starting  current  belonging  to  the  TE  13.1,1  with  3.2 
amps,  and  the  highest  starting  current  corresponding  to  the  TE14.1  at  15  amps. 
Although  the  clean  spectrum  would  appear  to  allow  q  1  modes,  ihe  higher 
starting  currents  eliminate  the  possibility.  Scaling  the  TE13.11  to  the  TEm.i.j 
yields  Lr  12  amps. 

The  misleading  aspect  of  Figures  2.2a  and  2.2b,  is  the  representation  of  each 
mode  as  a  delta  function  in  B.  In  reality,  each  mode  may  be  excited  over  a  range 
of  B  as  predicted  by  equation  (2.29)  prior  to  specifying  x  —  rinjn.  There  are 
actually  large  regions  of  mode  overlap.  The  B  value  in  Figures  2.2a  and  2.2b 


at  -^min  is  representative  of  where  the  mode  should  be  found  experimentally,  but 
this  is  not  precisely  accurate.  The  xmin  does  not  coincide  with  optimum  efficiency 
which  will  be  discussed  in  Section  2.3. 

The  final  consideration  of  the  linear  theory  will  provide  a  more  reasonable 
view  of  the  mode  spectrum.  This  is  a  computer  code  developed  by  Kreischer  [47] 
which  provides  plots  of  starting  current  values  as  a  function  of  the  magnetic  field. 
The  program  is  based  on  the  formulation  outlined  in  this  section,  restricted  to 
the  thin  annular  electron  beam  with  no  velocity  spread.  An  arbitrary  longitudinal 
profile  may  be  chosen  and  the  code  allows  for  higher  harmonics.  Equation  (2.21) 
is  evaluated  over  the  desired  range  of  x  after  f(z)  and  n  have  been  specified. 

Figures  2.3a  and  2.3b  are  the  linear  code  computer  plots  of  the  second  har¬ 
monic  0.3.1  and  whispering  gallery  cavity  spectra  over  the  4.0  to  6.0  tesla  range. 
The  vertical  scale  is  now  linear  from  0  to  10  amps.  The  Gaussian  field  profile  was 
chosen  and  n  —  2.  To  satisfy  the  computer  program  input  requirements  in  con¬ 
sidering  all  possible  modes,  some  averaging  of  several  equation  parameters  such 
as  Q., and  was  required.  This  procedure  and  the  approximations  made 
in  the  derivation  of  the  Gaponov  equation  account  for  the  minor  1  T  differences 
between  Figures  2.2  and  2.3. 

The  broad  range1  of  mode  excitation  is  now  evidenl  for  the  0.3.1  cavity  in 
figure  2.3a.  Again  space  on  the  plot  did  not  allow  naming  all  of  the  modes,  but 
several  of  the  more  isolaied  and  lower  starting  current  modes  are  labelled.  As  an 
example  of  the  discussion  concerning  the  domination  of  one  rotating  mode  over 
the  opposite  rotation,  both  the  TE4  4  )  and  TE4,4<1(..)  modes  are  shown  at 


Linear  Code  Plot  of  the  Whispering  Gallery  Cavity  Mode  Spectrum 


about  4.35  tesla.  The  TE4i4ii(_)  lies  within  the  TE4i4i1(4)  excitation  range  at  a 
higher  JST,  and  should  not  be  experimentally  evident.  This  consideration  is  ap¬ 
plicable  to  a  limited  extent  between  adjacent  modes,  although  mode  enhancement 
or  suppression  by  adjacent  modes  was  not  considered  in  the  single  mode  analysis 
presented  here.  Likely  examples  of  this  phenomenon  are  the  TE6,3,i(_)  and  the 
TE()  3  j(T)  which  are  the  unlabelled  lines  lying  within  the  TE9)2ii(-)  excitation 
zone  at  about  4.17  tesla.  As  a  first  approximation,  the  TEg.3,1  mode  would  not 
be  expected  to  appear,  though  it  is  shown  above  threshold  on  Figure  2.2a. 

The  whispering  gallery  plot  shown  in  Figure  2.3b  again  shows  the  reduced 
mode  excitation  spectrum.  The  seven  modes  predicted  to  be  above  threshold  by 
the  Gaponov  calculations  are  shown  as  well  as  the  TEr^.if-j  )•  The  discrepancies 
between  the  starting  current  values  in  the  figure  and  the  values  in  Table  2.2b  are 
accounted  for  by  the  parameter  averaging  procedure  which  was  used  to  produce 
computer  plots  of  7>x  for  all  the  modes  on  the  same  graph. 

The  mode  excitation  calculations  presented  are  the  most  useful  aspects  of 
the  linear  theory.  However,  theoretical  predictions  of  efficiency  and  output  power 
require  the  nonlinear  theory  applicable  beyond  the  t/  —  S 1  regime. 

2.3  The  Nonlinear  Theory  of  the  Gyrotron 

Gyrotron  operation  at  optimum  detuning  and  current  is  generally  well  beyond 
the  limits  of  applicability  for  the  linear  analysis.  Realistic  estimates  of  output 
power  depend  on  a  nonlinear  theory. 

The  slow  time  scale  analysis  is  based  on  the  work  of  Soviet  theorists.  The 
general  derivation  is  outlined  by  Flyagin  in  48  although  the  formulation  appeared 
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earlier  in  [7j.  The  theory  leads  to  two  dimensionless  coupled  differential  equations 
relating  the  particle  dynamics  and  the  electromagnetic  field.  The  equation  de¬ 
scribing  the  electron  motion  in  a  combined  RF  and  static  magnetic  field  in  the 
weakly  relativistic  limit  is  credited  to  Yulpatov  in  several  Soviet  papers  includ¬ 
ing  |49|.  B.  Danly  of  the  MIT  Plasma  Fusion  Center  rederived  the  equations  in 
terms  of  the  Soviet  parameter  definitions.  His  development  is  reconstructed  in  the 
following  subsection. 

The  difficulties  inherent  in  the  nonlinear  analysis  necessitate  a  single  particle 
approach.  This  is  equivalent  to  a  delta  function  for  the  electron  velocity  distribu¬ 
tion  function  as  discussed  in  the  linear  theory.  Other  restrictions  include  a  weakly 
relativistic  beam,  again  only  retaining  the  ^  in  the  denominator  of  the  cyclotron 
frequency.  Additionally,  space  charge  effects  will  be  neglected. 

—  Derivation  of  the  Slow  Time  Scale  Equations  — 

The  coordinate  system  remains  the  same  as  for  the  CAVRF  description,  with 
the  2-axis  aligned  with  the  static  magnetic  field  and  the  gvrocenter  as  the  origin 
in  the  x-y  plane.  The  Lorentz  force  equation  for  an  electron  in  the  RF  and  static 
magnetic  field  is 

~  -  -e[f  4(r  «  £)]  .  (2.34) 

This  equal  ion  is  separated  int  o  perpendicular  and  parallel  components  with  respect 
to  the  2-axis.  Neglecting  the  RF  magnetic  field,  the  B  in  equation  (2.34)  represents 
the  static  magnetic  field  only.  The  electric  field  component  £z  =  0  for  TE  modes, 


so  the  parallel  component  equation  yields  p n  =  constant.  The  equation  describing 
motion  in  the  transverse  plane  is 

-r  uc  (px  x  z)  =  -ei  (2.35) 

where  p±  is  the  nonrelativistic  perpendicular  momentum  mvx-  The  position  and 
momentum  vectors  may  be  written  in  terms  of  complex  variables  R  and  P  defined 

iP 

p„  --  P  -  p_e'^  and  r  -  R  -  *’•' 2)  =>  R  = - - —  .  (2.36) 

TTt  gOJ  c 

The  7T  2  phase  difference  between  R  and  P  is  representative  of  the  circular  motion 
about  the  gyrocenter.  The  perpendicular  equation  of  motion  (2.35)  is  now  written 

~  -  ijjcP  =  -ti  .  (2.37) 

at 

The  fields  must  be  converted  to  a  complex  notation.  With  the  condition  of  a 
slowly  varying  radius  along  the  z-axis  condition,  the  fields  are  separated  into 
longitudinal  and  transverse  components  by  if(i,t/,z,<)  -  Re  f(z)E?(x.y)e'ui1  . 
For  TK  modes,  the  transverse  field  components  ( E. )  are  derivable  from  the  scalar 
fund  ion  'F.  . 

Ef  z  ■  V  'F,  where  'F,  satisfies  V2  vF..-  —  A-2  ty..  -  0  (2.38) 

The  function  vF..  may  be  expanded  in  a  Fourier  series  in  terms  of  an  angular 
variable  6  about  the  axis  of  gyration. 


Using  Flyagin’s  notation,  the  part  of  the  field  moving  in  synchronism  with  the 
electron  and  therefore  interacting  the  most  efficiently  is  Esynch  which  rotates  at 
the  rtth  harmonic  of  uc. 


Esynch  Z  X 

Here  the  is  given  by  which  may  be  expressed  as 


(2.40) 


=  Anrne~'nfi  =  An(R')n  . 


(2.41) 
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The  R  in  the  equation  above  is  the  complex  position  coordinate  defined  in  (2.36). 
The  field  phase  8  measured  from  the  z-axis  may  coincide  with  an  arbitrary  initial 
electron  phase  p  n  2  because  there  will  be  an  ensemble  average  over  p  in  the 
efficiency  determination.  The  complex  form  of  the  equation  (2.38)  combined  with 
equation  (2.41)  yields 

Esynch  -  -  -f,l„V(f?T  •  (2.42) 

The  goal  is  to  find  an  applicable  complex  form  of  £  (x.y,  z,t)  written  as  f(z)E'stlMi. 
With  f{z)  real.  E'  must  be  chosen  to  satisfy  the  following  two  conditions  which 
lend*  to  1  he  stated  result . 

£'•>-  f{z)Rc(E[.e'^)  ) 

f  >  E[  -  Esynch -X  =  -iA„n(R-)n~1  (2.43) 

S'  -  y  f(z)hn(Ey-'<)  ) 

Writing  the  position  variable  R  in  terms  of  P  by  the  relation  (2.36)  and  substi¬ 
tuting  in  for  e.£  on  the  right  side  of  equation  (2.37)  produces  the  equation 


(meu:c)n 


at 
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(2.44) 


«-*-!*  •  »  ■  •  «  »  *  •  •  »  _  •  '  *  »  '  M  ***«'»*(*»’  •*«*»*»*  k^1  ■*'  ^  k^ 


%  v  o*  >"  ^  y\V»v vy^.  ^ "*.' v : 


The  complex  momentum  P  may  be  converted  to  the  slow  time  scale  momen¬ 


tum  p  by  the  relation 


P  —  pe 


+  i(u/n)t 


(2.45) 


Combining  this  with  Flyagin’s  definition  of  f 


einAnn 
(meu,'c) n 


(2.46) 


yields  the  slow  time  scale  equation 


(n  '  ^C)P  ^  *  ' 


(2.47) 


The  equation  is  written  in  dimensionless  form  by  converting  to  the  dependent 
variable  q  and  the  independent  variable  c  defined  by 


and  f  = 


2  v: 


(2.48) 


The  dimensionless  q  represents  the  momentum  normalized  to  its  value  at  the  start 
of  the  interaction  region.  All  of  the  “(<■)”  superscripts  indicate  the  parameter 
value  prior  to  the  interaction.  This  is  a  necessary  distinction  since  is  changing 
throughout  the  cavity.  The  normalized  longitudinal  coordinate  c  is  related  to  // 


defined  in  equation  (2.27)  by  <>*  =  ( zjLeff)p . 

Kqtialion  (2.47)  in  terms  of  these  dimensionless  variables  is 


3u')2u,!.p) 


a  -  -•) 


q  -  i(q' 


(2.49) 


where 


i^p{,,)  "■  2 


it-'Z  i.) 


(2.50) 


‘  v.vl" w.  v  .v.v  :  /.v;  c;  „y 


•  •  •  •  Ov'  .  *  •  - 

»  -  ■  m  »  .  r  ^  -  y  -  »  ■ 


The  bracketed  term  on  the  left  side  of  equation  (2.49)  may  be  simplified  by  the 
following  approximations  and  the  definition  of  the  Soviet  detuning  parameter  A. 


Since  p\\  is  constant,  'j^/3 ^  =  10\\.  The  weakly  relativistic  approximation  com¬ 
bined  with  u  ss  nuc  allows 

r 

2 


fi[o)2Jeo) 


(;-»*) 


A+  \q\z  -  1  . 


(2.51) 


The  A  is  a  measure  of  the  off  resonance  condition  and  is  defined  explicitly  as 


A  - 


u.'  —  nu: 


(o) 


3 


(o)2 


A  - 


2x 


(2.52) 


The  statement  after  the  arrow  relates  A  to  the  detuning  parameter  x  and  the 
normalized  cavity  length  p  introduced  previously. 

The  derivation  is  now  complete  to  present  the  Yulpatov  equation  in  the  form 
appearing  in  [48], 


dq 


:  +  *(A  +  l9i  -  i )q  =  *(<?")"  ? 


(2.53) 


This  equation  describes  the  electron  dynamics  during  the  interaction.  To  describe 
the  fields,  Maxwell's  equations  may  be  used  to  arrive  at 


d2f  ,  ,2 


.V  idX 


k  f  -  -  /  dxdyJ^.  ■  E: 


(2.54) 


where 


2tt 


A',,.  -  ~  I  dxdy"$> .J2  and  Ju.  —  —  J  J(t)e  luJld(^t) 


(2.55) 


Assuming  an  infinitely  thin  beam  at  r  =  re,  the  transverse  current  density  may 
be  written  in  terms  of  the  charge  density  p  =  -1  jv,\. 


J±(t)  =  P  —  f>(r  -  re) 
rn, 


(2.5(») 


57 


a!* 


This  combined  with  Es  —  E's  given  in  equation  (2.43)  yields 


”  jdxdyJl  ■  J 5;  =  X  •  (2.57) 

Here  the  represents  the  integral  of  the  bracketed  quantity  with  respect 

to  ujt  over  a  period,  divided  by  2n.  The  substitution  t  — ►  t0,  where  t0  represents 
the  time  that  the  interaction  begins,  is  permissible  because  the  r  determining  p 
is  a  constant.  The  \  in  equation  (2.57)  is  defined  as 


\  = 


Ss[im 


(2.58) 


Equation  (2.54)  now  becomes 


k?{f  =  x(Pne 


—  I  u  t  \ 

■■  wt.. 


\  p 


(2.59) 


the  last  step  a  result  of  converting  to  the  slow  lime  scale  defined  by  equation  (2.45). 
All  that  remains  is  changing  to  the  normalized  variables  q  and  s*.  Omitting  the 
algebra,  the  resulting  equation  is 


(is1 


°‘f 


(2.00) 


lleie  the  dimensionless  current  i  and  c  are  normalization  (actors. 


2  k -v  Heir  . . 

a  and  I  - - -  -  v~  .4,.  *  r("  *1,;  (2.01 

3{"V1 


The  r[i,)  in  the  definition  of  I  is  the  siow  time  scale  equivalent  of  R  defined  in 
equation  (2.36).  The  primary  results  of  the  slow  time  scale  derivation  are  now 
restated. 


dq 


i  ?( A  4  | q  -  -  1)7  -  t(<7')"~  1  J 


(2.02) 
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V. 


(2.63) 


='<«■> 

These  two  equations  are  coupled  by  the  /  which  is  in  the  definition  of  7  in 
the  first  equation  and  the  explicit  appearance  of  q  in  the  second.  Equation  (2.63) 
is  the  generalization  with  a  source  term  of  expression  (2.8)  solved  for  by  CAVRF. 
Since  q{ f)  represents  the  normalized  momentum  of  the  electron  as  it  traverses  the 
resonator.  q($on t)'2  represents  the  energy  lost  by  the  electron  at  the  completion 
of  the  interaction.  An  ensemble  average  over  the  input  phases  of  the  electrons  will 
produce  the  perpendicular  efficiency  by  the  relation 

Vi  -  1  -  J  dpo\q{<* ut);2  (2.64) 

0 

where  the  integration  variable  <p 0  represents  the  input  phase  in  the  boundary 
condition  q(°)  =  et<p".  The  electronic  efficiency  qee  and  the  output  efficiency  r)  are 
then  calculated  by  equations  (2.5)  and  (2.6). 

Application  of  the  Slow  Time  Scale  Theory  — 

The  intractable  nature  of  coupled  differential  equations  requires  a  numeri¬ 
cal  analysis,  and  most  theories  make  further  approximations  to  reduce  equations 
(2.62)  and  (2.63).  A  particularly  useful  approach  was  formulated  by  G.  S.  Nusi- 
novich  and  R.  E.  Erin  8  .  Their  simplification  was  based  on  the  premise  that 
for  high  Q  cavities,  the  longitudinal  profile  function  f{z)  is  fixed  by  the  cavity 
profile.  For  most  applications,  the  Gaussian  function  specified  by  equation  (2.12) 
will  provide  the  approximation  for  f(z).  In  terms  of  the  normalized  parameters 

m  =  e  .  (2.65) 
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This  approximation,  discussed  extensively  in  deriving  the  Gaponov  equation  in 
Section  2.2,  eliminates  the  consideration  of  equation  (2.63)  altogether.  The  result 
is  a  single,  first  order  differential  equation  which  may  be  integrated  numerically. 

The  equation  treated  by  Nusinovich  and  Erm  was  of  the  same  form  as  (2.62), 
but  differed  slightly  in  the  definitions.  The  exact  equation  listed  in  their  1972 


paper  is 


-  -  (A  +  .a|2  -  l)a  =  ian  1  f($)F  . 


dc  n 


(2.66) 


The  definitions  of  c.  //  and  A  remain  the  same,  but  the  F  which  appears  in  the 


same  position  as  Flyagin's  J  from  (2.50).  is  defined  by 


n-4  (  \ 

\2 ) 


J m~n  (k jl  f  e) 


(2.67) 


The  field  equation  equivalent  to  (2.63)  is  not  listed  in  their  paper  because  they 
have  assumed  a  Gaussian  for  /(c),  but  the  normalization  leads  to  a  definition  of 
their  dimensionless  current  /,,  different  from  the  1  appearing  in  Flyagin's  analysis. 


I  0.2-1  I  (Q  ■  10  ')  7t 


2(3  -it l 


(¥)' 


(2.68) 


I  lie  (i  i«-  the  same  beam  interaction  term  defined  b\  equation  (2.28)  which  appears 
in  t  he  Gaponov  linear  theory  equation.  The  expression  for  /,.  is  dimensionless  for 
a  beam  current  /  in  amps. 

The  results  of  the  work  by  Nusinovich  and  Erm  were  plots  of  r)±  and  optimum 
detuning  A.,,,/  over  ranges  of  the  parameters  F,  p  and  I0.  This  was  accomplished 
by  specifying  a  value  of  F.  //  and  A.  then  integrating  equation  (2.66)  over  the  range 
of  c  from  to  -f  -F’fi.  The  parameters  were  held  fixed  while  the  integration 


was  performed  32  times  with  separate  input  phase  <p0  values.  This  allowed  a 
numerical  integration  of  equation  (2.64)  and  calculation  of  r?_.  The  procedure 
was  repeated  varying  only  the  A  in  equation  (2.66)  until  the  value  of  A  which 
produced  the  highest  efficiency  (A opt)  was  located.  By  running  the  algorithm  for 
many  values  of  F  and  //,  isoefficiency  curves  and  the  values  of  Aopt  were  plotted 
against  F  and  //.  At  each  value  of  F,  n,  and  corresponding  to  A opt,  the  value 
of  1K,  was  computed  by  the  equation 

1  f,  =-  FV  V(3-n)  ,  (2.69) 

which  is  derived  by  equating  the  RF  power  gain  and  the  power  given  up  by  the 
beam.  The  generated  graphs  of  ?/_  and  Ac,;)(  against  specific  values  of  Ic  and  fi 
are  extremely  useful  in  predicating  gyrotron  efficiencies  at  particular  operating 
conditions. 

Figure  2.4a  is  the  plot  of  r/_  and  A opt  as  a  function  of  I,,  and  p  which  resulted 
from  the  integrations  at  the  harmonic  value  of  n  —  2.  Figure  2.4b  is  the  same 
type  of  graph  but  differs  by  the  integration  length  used.  A  greater  extent  of 
the  Gaussian  profile  function  is  included  on  the  cavity  input  side  in  the  data 
represented  by  Figure  2.4b,  by  integrating  <r  from  —  \/3p  to  -t  -2"p.  Both  results 
are  credited  to  Xusinovich  and  Erm,  Figure  2.4a  is  from  their  1972  paper  [8j,  and 
Figure  2.4b  was  included  in  the  publication  by  Gaponov  et  al.  45:.  Both  figures 
have  linear  scales  on  the  horizontal  axis  for  //  and  log  scales  on  the  vertical  axis 
for  /,,.  but  caution  should  be  used  in  comparing  values  directly  because  the  relative 


scales  are  not  the  same. 


The  isoefficiency  curves  are  labelled  by  their  r) ±  values.  The  dashed  lines 
represent  the  A opt  value  which  yielded  the  given  t)x.  The  highest  efficiency  pre¬ 
dicted  at  the  second  harmonic  by  Figure  2.4a  is  72%  corresponding  to  lc  «  20  and 
p  a:  18.  The  highest  r]±  shown  in  Figure  2.4b  is  76%,  but  it  occurs  at  10  as  200 
and  p  a;  50.  The  higher  efficiency  in  this  region  by  the  second  graph  is  a  result  of 
the  longer  field  interaction  on  the  input  end  prebunching  the  electron  beam  and 
enhancing  efficiency.  However,  in  gyrotron  design  considerations,  the  long  cavity 
lengths  described  by  this  large  value  of  p  are  unattractive  because  Qo  ~  p2,  and 
the  higher  overall  Q  values  reduce  the  output  efficiency  rj  by  the  relationship  (2.5). 

In  practice,  the  specific  operating  conditions  and  cavity  length  determine  I0 
and  p  which  are  then  used  to  calculate  r)±  and  A0;,<  from  the  figures.  In  the 
present  study,  efficiency  calculations  were  made  for  all  of  the  second  harmonic 
modes  predicted  to  be  above  threshold  by  the  linear  analysis.  The  results  of  those 
calculations  are  presented  in  Tables  2.3a  2.4b.  The  shorter  integration  length  in 

Figure  2.3a  generated  the  values  presented  in  the  first  set  of  tables.  Table  2.3a  is 
for  the  0.3.1  cavity  and  Table  2.3b  is  the  whispering  gallery  cavity  data.  The  longer 
integration  length  represented  in  Figure  2.3b  was  used  to  compute  efficiencies  in 
Tables  2.4a  and  2.4b. 

The  values  of  p  ranged  from  15.8  to  21.7  for  the  modes  evaluated  in  the  0,3,1 
cavity.  The  whispering  gallery  modes  with  a  shorter  Lci\.  produced  p  values  from 
10.5  to  13.6.  The  parameter  was  calculated  by  equation  (2.68)  for  a  beam 
current  of  5  amps.  As  shown  in  the  tables,  the  resulting  values  of  la  ranged 
from  0.4  to  3.2  in  the  0.3.1  cavity  and  from  0.52  to  0.81  in  the  whispering  gallery 
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Table  2.3a  —  0,3,1  Cavity  uj  rs  2u>c  Efficiency  Calculations 
Nusinovich  and  Erm  —  1972 

/o  /  Q 

Integration  Length  of  — f  M  <  ?  <  ty/J 


MODE 

h 

F 

&  opt 

n±  (%) 

n(%) 

TE1,5,1<4) 

15.79 

1.3 

0.042 

0.31 

33 

20 

TE63i,(4  ) 

16.13 

1.6 

0.047 

0.33 

37 

22 

TE9,2,1(-1 

16.15 

2.5 

0.066 

0.37 

46 

26 

te4  4  ,  i  h  j 

16.72 

0.81 

0.026 

0.18 

23 

14 

TE2.5.1(, ) 

17.04 

1.6 

0.045 

0.32 

37 

22 

TE10  2.i( -  ) 

17.13 

2.5 

0.063 

0.37 

46 

25 

TEo.6,1 

17.15 

0.90 

0.027 

0.19 

24 

14 

TE7  3  1(4  ] 

17.19 

1.4 

0.041 

0.30 

35 

20 

TE5  4  1(4  ) 

17.85 

1.5 

0.041 

0.31 

36 

21 

TEj  ]  .2,3  ( -  1 

18.10 

2.4 

0.057 

0.35 

44 

23 

TEg .3,1(4 ) 

18.24 

1.2 

0.034 

0.30 

33 

18 

TE3.6.1(+  ) 

18.25 

1.4 

0.039 

0.30 

36 

21 

TE6  4  1(  4 ) 

18.97 

2.2 

0.051 

0.34 

42 

24 

TE12.21(-  ) 

19.06 

2.1 

0.049 

0.33 

42 

21 

TE9  3  1(4.  ) 

19.28 

0.99 

0.027 

0.18 

28 

15 

TE4.5i,(4  ) 

19.44 

0.96 

0.024 

0.25 

24 

14 

TE2.6  1(- ) 

19.71 

0.77 

0.021 

0.24 

24 

13 

TE0  6  1 

19.80 

0.40 

0.011 

0.15 

18 

9.9 

TE„.2.1(.  ) 

20.02 

1.7 

0.040 

0.31 

38 

19 

TE7  4  1  (  —  1 

20.07 

2.8 

0.056 

0.35 

46 

25 

TEio.3.1,-  ) 

20.31 

0.79 

0.024 

0.20 

30 

16 

TE5.6.K- ) 

20.60 

1.7 

0.039 

0.29 

38 

21 

TE3  6  ]  (  | 

20.94 

1.8 

0.039 

0.29 

38 

21 

TE]4  2,|(-  ) 

20.98 

1.3 

0.032 

0.23 

34 

16 

TE],7,h_  j 

21.10 

1.7 

0.038 

0.29 

37 

20 

TEg  4  1(4  | 

21.16 

3.2 

0.053 

0.35 

48 

25 

TE11.311  ~  j 

21.32 

17 

0.037 

0.29 

38 

19 

TE6  6,1(-) 

21.66 

1.8 

0.038 

0.31 

38 

21 

19 

21 
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Table  2.4a  —  0,3,1  Cavity  u;  ss  2wc  Efficiency  Calculations 
Nusinovicb  and  Erm  —  1975 

Integration  Length  of  -s/Sfi  <  f  <  + 


MODE 

h 

F 

^ Opt 

>?,  (%) 

n(%) 

TEi,5.1(t  ) 

15.79 

1.3 

0.041 

0.30 

32 

19 

TE6,3,1(-*-) 

16.13 

1.6 

0.046 

0.31 

35 

21 

TE9i21(_  ) 

16.15 

2.5 

0.064 

0.35 

43 

24 

TE4  4.J(.  | 

16.72 

0.81 

0.026 

0.19 

23 

14 

TE2.5.1(  -  ) 

17.04 

1.6 

0.044 

0.30 

36 

21 

TE102.I(  | 

17.13 

2.5 

0.061 

0.35 

43 

23 

TEo.5.1 

17.15 

0.90 

0.027 

0.23 

24 

14 

TEr.s.„, ) 

17.19 

1.4 

0.040 

0.29 

34 

19 

TEs.41H  ) 

17.85 

1.5 

0.040 

0.29 

34 

20 

TEj  j  2 ,u-  ) 

18.10 

2.4 

0.055 

0.34 

41 

22 

TE8,3,1(-  ) 

18.24 

1.2 

0.033 

0.25 

31 

17 

TE3.6.1(+) 

18.25 

1.4 

0.038 

0.28 

34 

20 

te6i4i1M 

18.97 

2.2 

0.049 

0.32 

40 

22 

TE12  21(.  ) 

19.06 

2.1 

0.047 

0.31 

39 

20 

TE9  3  jj-  ) 

19.28 

0.99 

0.028 

0.17 

30 

16 

TE4.6.„.  ) 

19.44 

0.96 

0.027 

0.29 

28 

16 

TE2  6  J(.  > 

19.71 

0.77 

0.023 

0.19 

28 

16 

TE0  6  i 

19.80 

0.40 

0.016 

0.13 

26 

14 

TE]32.1(  i 

20.02 

1.7 

0.039 

0.30 

36 

18 

te7  4  i i ..  j 

20.07 

2.8 

0.054 

0.33 

42 

23 

TEjo  3  1( .  j 

20.31 

0.79 

0.023 

0.18 

28 

15 

TEs.6.I(-  ) 

20.60 

1.7 

0.038 

0.29 

36 

19 

TE3i6,1(- ) 

20.94 

1.8 

0.038 

0.28 

36 

20 

te14  21(  j 

20.98 

1.3 

0.031 

0.21 

32 

15 

TEj ,7,i (-  ) 

21.10 

1.7 

0.036 

0.28 

35 

19 

TE8|4,1(,) 

21.16 

3.2 

0.055 

0.33 

43 

23 

TEjj  3 

21.32 

1.7 

0.036 

0.27 

35 

18 

te,,m-) 

21.66 

1.8 

0.036 

0.29 

35 

19 
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Table  2.4b  —  Whispering  Gallery  Cavity  w  ss  2wc  Efficiency  Calculations 

Nusinovich  and  Erm  •—  1975 

Integration  Length  of  \  3/j  <  s-  ^  ^ 


MODE 

n 

/,> 

F 

A,,, 

^  (%) 

r>(%) 

TE9,i,i(-) 

10.56 

0.81 

0.026 

0.28 

9.1 

5.5 

TEio,i,i(-  ] 

11.24 

0.81 

0.028 

0.27 

12 

7.2 

TEn  j  i,  , 

11.92 

0.80 

0.028 

0.26 

14 

8.1 

TEj2.i.](  ) 

12.61 

0.77 

0.027 

0.25 

15 

8.3 

te82.J(  .  j 

12.77 

0.52 

0.010 

0.20 

3.1 

2.0 

TE131  J(  1 

13.31 

0.73 

0.025 

0.25 

15 

8.2 

^E9.2.J  (  -  ) 

13.55 

0.53 

0.013 

0.20 

6.0 

3.8 

cavity.  From  the  isoefficiency  plots,  it  is  evident  that  the  optimum  efficiency  at 
the  second  harmonic  in  the  10  <  <  20  range  occurs  at  values  of  10  ss  20.  Since 

1,  ~  7,  the  optimum  efficiency  would  require  an  order  of  magnitude  increase  in 
beam  current. 

These  low  values  of  10  which  model  the  experimental  investigation  require  a 
large  degree  of  interpolation  in  using  the  graphs.  The  problem  is  most  obvious  for 
the  ],,  1  situations,  which  would  plot  below  the  horizontal  axis  in  Figures  2.4a 

and  2.4b.  The  solution  was  a  computer  code  developed  by  B.  Danly  at  the  MIT 
Plasma  Fusion  ('enter  which  solved  the  differential  equation  (2.66)  in  terms  of  the 
same  parameters  defined  in  (2.67)  (2.69).  There  are  two  versions  of  the  code 

to  allow  investigation  of  both  integration  lengths  discussed  previously,  and  32 
steps  are  included  in  the  integration  over  input  phases.  The  program  reproduces 
the  Nusinovich  and  Errri  results,  but  is  more  flexible  by  allowing  investigation 
over  selected  ranges  of  7\  p  and  with  results  available  in  either  tabulated  or 
graphic  form.  This  code  produced  all  of  the  values  of  r/  ^  listed  in  Tables  2.3a  - 
2.1b.  The  last  column  in  each  table  lists  the  output  efficiency  r/.  related  to  i]± 
by  equations  (2.5)  and  (2.6). 

The  largest  value  of  //  listed  in  Table  2.3a  for  the  0.3.]  cavity  is  the  TEbi4j(_) 
at  48'','  with  the  lowest  efficiency  belonging  to  the  TEo.0,1  with  r/  -  18*7.  In 
contrast,  the  whispering  gallery  cavity  results  for  the  short  integration  length 
show  2.4*7  •  7/  <  14*7.  This  difference  in  second  harmonic  mode  efficiencies 

would  be  expected  from  the  linear  results,  where  the  stronger  modes  in  the  0,3,1 
cavity  had  7  T  ^  0.5  amps  and  for  the  whispering  gallery  cavity  l.T  &  3.5  amps. 


Comparison  of  the  results  between  the  two  integration  lengths  reveals  that 
the  differences  in  r\ _l  are  generally  limited  to  less  than  3%.  Specific  trends  indicate 
that  the  isoefficiency  contours  are  more  closely  spaced  for  the  shorter  integration 
length  producing  a  greater  range  in  the  values  of  r)±.  The  two  lengths  produce 
the  same  results  at  ss  30%.  For  >  30%.  the  values  tabulated  in  Table  2.3a 
are  greater  than  those  in  2.4a  and  with  r? _  <  30%.  the  opposite  is  true.  The 
trend  is  supported  by  the  whispering  gallery  cavity  results  which  are  in  the  low 
efficiency  regime,  consequently  the  values  of  i]  in  Table  2.3b  are  low'er  than  those 
in  2.4b.  This  analysis  is  valid  only  over  the  ranges  of  J,  and  p  covered  by  the 
tables.  Figures  2.4a  and  2.4b  show  that  the  relative  contour  spacing  between  the 
short  and  long  integration  lengths  varies  between  different  regions  of  10  and  //. 

As  a  final  consideration  of  the  slow  time  scale  nonlinear  theory,  it  should  be 
mentioned  that  Fliflet  and  Read  } 50)  have  developed  a  self-consistent  numerical 
solution  to  equations  equivalent  in  content  to  (2.62)  and  (2.63).  By  solving  the 
coupled  equations  with  an  iterative  procedure  until  a  specified  tolerance  is  reached, 
ihe  resulting  solution  should  yield  increased  accuracy  over  the  simple  assumption 
of  a  fixed  functional  form  of  f{z).  The  code  was  applied  to  ihe  0,3.1  cavity 
TK]|,2.i(  )  mode  in  the  linear  limit,  but  a  large  discrepancy  between  the  results 
and  the  Gaponov  and  linear  code  predictions  cast  doubt  that  the  program  was 
executing  properly  at  the  second  harmonic. 

-  Fast  Time  Scale  Efficiency  Calculations  — 

In  the  derivation  of  equations  (2.62)  (2.63).  a  slow  time  scale  momentum 

variable  p  was  introduced  as  defined  in  (2.45).  This  variable  eliminated  the  change 


in  momentum  due  to  the  cyclotron  rotation.  Therefore  dp/dt  =  0  in  the  absence 


of  any  interaction  with  the  RF  field. 


A  computer  code  in  use  at  the  MIT  Plasma  Fusion  Center  integrates  the 


particle  equations  of  motion  in  terms  of  the  fast  time  scale  momentum  P  intro¬ 


duced  in  equation  (2.36).  Initially  developed  by  K.  R.  Chu  of  the  Naval  Research 


Laboratory,  the  program  also  differs  from  the  derivation  of  the  slow  time  scale 


equations  by  including  the  effects  of  the  RF  magnetic  field.  Because  the  code  in¬ 


tegrates  along  the  actual  particle  trajectories  with  the  full  interaction  of  the  high 


frequency  fields,  it  is  referred  to  as  a  panicle  integrator.  The  penalty  for  the  more 


complete  system  of  equations  is  the  increased  computer  processing  time  required. 


Chu  eliminated  the  need  for  the  field  equation  (2.63)  by  assuming  a  sinu¬ 


soidal  profile  for  f[z).  This  differs  from  the  approach  of  Nusinovich  and  Erm 


who  modeled  f{z)  with  a  Gaussian.  A  significant  modification  of  the  original  fast 


time  scale  code  was  developed  by  J.  Schutkeker  of  MIT.  The  new  version  accepts 


the  actual  cold  cavity  CAVRF  field  profile  as  a  tabulated  function  for  f{z).  The 


perturbing  effect  on  the  fields  by  the  electrons  is  therefore  not  included,  but  the 


CAVRF  profile  should  provide  a  more  accurate  representation  for  f(z)  then  either 


a  Gaussian  or  sinusoidal  function. 


The  particle  integrator  requires  as  input  a  normalized  version  of  the  CAVRF 


output,  the  detuning  parameter  x  defined  in  (2.20).  and  the  normalized  field 


strength  E.  In  terms  of  the  parameters  defined  in  the  Nusinovich  and  Erm  theory, 


FH.i'  "  /  2"  1  n! 


./»,!»(*'  r,  )  \  tin 


(2.70) 


v  V.v.v.^v- . .  C.  -v-v  .  -  . 


■* f 


Here  the  B  is  a  normalized  magnetic  field  given  by  B  =  ^nu;cr0/c.  The  F  is 
defined  in  equation  (2.67)  and  may  be  used  to  relate  E  to  the  field  amplitude  A, 
with  the  result  E  =  2A/n. 

In  execution,  the  program  separately  evaluates  each  input  value  of  x  and  E. 
The  code  uses  the  CAVRF  diffractive  Q  value  and  E  to  compute  an  output  power 
by  Pout  ~  E2/Q[).  The  input  power  required  to  produce  this  Poui  is  determined 
by  the  T)el:  value  from  the  integration  along  the  electron's  trajectory.  The  final  step 
simply  solves  for  the  beam  current  needed  to  produce  Pm  at  a  specified  voltage. 

The  correct  value  of  E  at  optimum  detuning  xopt  must  be  found  which  pro¬ 
duces  the  value  of  /  representing  experimental  operating  conditions.  The  xopt  is 
optimized  with  respect  to  r)c(  for  a  given  value  of  E.  A  reasonable  estimate  of  the 
correct  value  of  E  is  given  by  first  performing  the  calculations  in  the  Nusinovich 
and  Erm  slow  time  scale  theory,  then  using  the  resulting  F  value  in  equation  (2.70). 
In  a  similar  manner,  x(,vt  is  related  to  A bv 

*o,.t  -  -  '2r'  V  •  (2-71) 

These  estimates  of  E  and  x„vi  will  be  accurate  only  to  the  extent  that  the  calcu¬ 
lated  efficiencies  between  the  two  theories  agree. 

As  a  verification  of  the  accuracy  of  the  particle  integrator  results  at  the  sec¬ 
ond  harmonic,  the  code  was  run  in  the  linear  limit  for  the  whispering  gallery 
TK|2.i,i(  ,  mode.  The  linear  limit  was  confirmed  by  inputting  several  low  values 
of  E  and  checking  that  ?/  —  E2.  The  current  at  these  values  is  equal  to  /-T.  To 
compare  with  the  Gaponov  results,  a  Gaussian  function  was  used  for  f(z)  instead 
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of  the  CAVRF  profile  function.  The  only  difference  between  the  theories  under 
these  conditions  should  be  the  consideration  of  the  RF  magnetic  field  by  the  parti¬ 
cle  integrator.  The  agreement  was  very  good,  with  the  Gaponov  equation  yielding 
/.- T  -  3.3  amps  and  the  fast  time  scale  code  producing  an  answer  of  /ST  =  3.45 
amps.  Additionally,  the  Gaussian  function  for  f(z)  produced  an  overall  efficiency 
calculation  of  r/  =  7.1%.  at  5  amps  for  the  same  mode.  This  is  close  to  the  slow 
time  scale  result  of  rj  -  8.1(7  shown  in  Table  2.3b. 

Efficiencies  were  evaluated  for  ten  inodes  in  the  0.3.1  cavity  using  the  full 
capabilities  of  the  fast  time  scale  code  with  the  CAVRF  longitudinal  field  distri¬ 
bution  and  two  separate  integration  lengths.  The  values  presented  in  Table  2.5a 
resulted  from  using  an  integration  length  of  —^n  S  c  <  -  *,'•  similar  to  Nusi- 
novich  and  Erm  1972.  The  longer  integration  length  of  -2v'3p  <  ?  <  -r2v//3p 
produced  the  efficiencies  listed  in  Table  2.5b.  Besides  the  predicted  values  of  r)e ( 
and  the  corresponding  values  of  E  and  xovt  are  shown.  The  results  correspond 
to  l  ,  04.3  k\’  and  15  amps.  The  long  integration  length  efficiencies  were 

slightly  lower  than  the  short  length  results,  bul  the  number  of  modes  evaluated 
precludes  reaching  a  general  conclusion. 

Tabic'  2.0  compares  the  slow  time  scale  and  fast  time  scale  results  for  these 
ten  modes.  The  particle  integrator  consistently  calculated  lower  efficiencies  than 
the'  Nusinovich  and  Erm  theory  predictions.  The  relative  differences  between  the 
computed  values  varied  between  the  modes  selected.  The  TE]  1,2,1  ( -)•  TEi2,2.i ( ..  j 
and  TE|:»>2.i(  )  showed  the  best  agreement  and  the  largest  disagreements  were 
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Table  2.5b  -  0.3,1  Cavity  ~  2u/’,  Efficiency  Calculations 
Fast  Time  Scale  Particle  Integrator 

Integration  Length  of  2\'3/i  < '  <;  <  ->  2\/3 fj 


MODE 

E 

x»v< 

T)"  (c?) 

ri  {%) 

0.30 

2.3 

16 

14 

TE9  2  j(  _  ) 

0.44 

-2.8 

25 

21 

TE102.H  ^ 

0.49 

-2.7 

26 

20 

0.54 

-3.0 

25 

20 

TE12.3i„.| 

0.58 

-2.9 

25 

18 

te13,3J(  ) 

0.58 

-2.5 

21 

15 

te7  4  ](4 ) 

0.48 

-2.7 

18 

15 

TEb.b.k- ) 

0.53 

-2.6 

22 

18 

TE3_e.11  - ) 

0.54 

-  3.0 

21 

17 

TEfi  4.1  ( .  j 

0.59 

-2.9 

23 

18 

p 

IL 
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Table  2.6  —  0,3.1  Cavity  w  ss  2 wc  Efficiency  Calculations  Comparison 

EFFICIENCY  -  t)  (%) 

NVSINOYICH  and  ERM  PARTICLE  INTEGRATOR 


MODE 

SHORT 

LONG 

SHORT 

LONG 

TE,.5  ,(  +  ) 

20 

19 

16 

14 

TEg  2  ](_) 

26 

24 

22 

21 

TEjo.a.n- ) 

25 

23 

21 

20 

TEii?2iJ(_  ; 

23 

22 

20 

20 

TE12i3.n(_, 

21 

20 

18 

18 

TEi3.2,l(-  ) 

19 

18 

17 

15 

TE7.4.K-  ) 

25 

23 

15 

15 

TEb.s.k  ) 

21 

19 

18 

18 

TEjc k  i 

21 

20 

18 

17 

TE8  4  !(..  I 

25 

23 

18 

18 

with  the  TEi  51(+)  and  the  TEy^,^)  modes.  The  particle  integrator  calcula¬ 
tions  for  the  whispering  gallery  cavity  continued  the  trend  of  lower  efficiencies  by 
predicting  all  of  the  modes  to  be  below  threshold.  Because  the  fast  time  scale 
program  agreed  reasonably  well  with  the  slow  time  scale  theories  when  input  with 
a  Gaussian,  the  major  source  of  the  discrepancy  was  believed  to  be  the  difference 
in  field  distribution  between  the  Gaussian  and  the  CAVRF  profile. 

2.4  vc\.  -  2-t'  Efficiency  Scaling 

The  slow  time  scale  nonlinear  theory  with  the  Gaussian  approximation  for 
j(~)  suggests  a  method  of  estimating  efficiencies  at  the  second  harmonic  from 
the  operating  parameters  at  the  fundamental.  The  simplicity  of  the  efficiency 
determinations  using  the  Nusinovich  and  Erm  plots  (Figures  2.4a  and  2.4b)  is 
the  result  of  having  reduced  r/j  to  a  function  of  the  two  parameters  10  and  fi¬ 
fty  deriving  a  scaling  law  between  /,,  at  the  fundamental  ij)  and  the  second 

harmonic  (/<,<„  -2>)  as  w°fi  a*  deriving  the  relationship  between  p(„  j)  and  //(„_ 2), 
the  graphs  may  be  used  to  map  regions  of  interest  at  the  fundamental  to  the 
second  harmonic.  This  could  provide  a  reasonable  estimation  of  second  harmonic 
efficiencies  Irom  a  device  designed  for  operation  at  the  fundamental.  The  following 
analysis  was  originally  developed  bv  R.  J.  Ternkin  at  MIT. 

The  parameters  necessary  for  the  scaling  aproxirnation  are  those  in  equa¬ 
tion  (2.27)  defining  //  and  equation  (2.fiK)  for  /,..  The  derivation  will  be  kept  as 
general  as  possible,  though  some  assumptions  will  be  made  concerning  the  oper¬ 
ating  mode  and  beam  placement.  All  of  the  experimentally  controlled  variables 


will  be  held  constant. 


Using  the  approximation  uj  «  nwc  ~  nB,  a  constant  magnetic  field  value 


yields  W(n=2)  ~  2w(n=1)  — *  A(n=2)  %  0.5A(n=j).  Neglecting  minor  differences 


between  the  effective  cavity  lengths  results  in  the  expression 


(*22)  *2(—) 

V  ^  /  (n  =  2)  V  ^  /  (n=l) 


(2.72) 


This  relation  and  cancelling  the  common  factors  in  equation  (2.27)  yields  the 
relationship  between  p  at  the  fundamental  and  the  second  harmonic. 


^ < tj — 2 )  ~  2p(T1=n 


(2.73) 


To  determine  the  scaling  law  for  /,,,  the  ratio  of  l0(n-~2)  A>(n=i)  's  formed  us¬ 
ing  equation  (2.68).  By  eliminating  the  constants  and  using  the  substitution  (2.72), 


the  following  expression  is  obtained. 


aril  sz  2.  ( il 
1)  ~  n2  \ 


)  (zT")  (2.74) 

/  V  ^eff  /  (n=i)  V(r!  =  l)  ^(n-i) 


To  proceed  further,  an  estimation  must  be  made  for  the  ratios  of  Q  and  G  between 


larmomcs. 


for  the  Q  scaling,  it  will  be  assumed  that  the  beam  placement  allows  ignoring 
the  whispering  gallery  modes.  From  previous  descriptions,  the  0.3.1  cavity  would 
meet  such  a  condition,  but  obviously  not  the  whispering  gallery  cavity.  Neglecting 
the  whispering  gallery  modes  permits  the  approximation  Q  ~  Qo .  The  expression 


for  Qi>  is  given  by 


where  R\  and  i?2  are  the  reflection  coefficients  for  the  ends  of  the  cavity.  The 
I-iflj  s  term  generally  be  slightly  larger  for  the  second  harmonic,  but  disre¬ 
garding  this  difference,  Q  p  ~  (^tt)2  yielding 


Q(r>=  1) 


(2.76) 


The  G  factor  defined  in  equation  (2.28)  is  obviously  mode  specific.  By  assum¬ 
ing  that  the  beam  will  be  in  a  near  optimum  location  to  excile  the  second  harmonic 
mode,  the  following  coarse  simplification  can  be  made  for  both  the  n  ~  1  and  n  -  2 
modes. 

Jl±n(kire)  _ 

^ml^mr) 

This  reduces  the  G  factor  scaling  to 


G(n  =  2) 
G(n  =  l) 


2 

mf>(n—  1 ) 


-  m 


2 

(n=l) 


"m,.(n  =  2)  ~  m(U2) 


(2.78) 


Having  already  disregarded  whispering  gallery  modes,  now  neglecting  other  ex¬ 
treme  "surface"  modes  allows  the  generalization  u~nj.  m-.  The  resulting  sim- 
plification  is  given  by 


(n  -  2) 

V2  ,  ,, 

k2 , 

A? 

mp(n~  1) 

~  L’h' 

(n=  1  ) 

^ m]i(n=  2) 

K±(n=2) 

~  A 2 

O'  - 

(2.79) 


This  scaling  factor  was  checked  against  actual  G  value  ratios  computed  using  the 
wall  and  beam  radius  values  applicable  to  the  0.3,1  cavity.  Modes  were  chosen 
such  that  =  2)  ==  („-]).  The  results  produced  a  range  of  G ..-j,  G^„  .  n 

ratios  from  0.12  to  0.46.  It  appears  that  the  scaling  is  accurate  to  within  a  factor 
of  two  with  the  included  mode  restrictions  and  beam  placement  assumptions. 


Expressions  (2.76)  and  (2.79)  for  the  Q  and  G  ratios  cancel  each  other  when 
substituted  back  into  equation  (2.74).  The  remaining  factors  are  easily  written  in 
terms  of  /i(n=J). 

I/4WM2  Ml 

X2\3±)  Ueff  /(„  =  ,)  M(rt- 1} 

The  final  expression  for  IQ{n= 2)  in  terms  of  7,,  at  the  fundamental  is 

2:7  2 

^o(n  —  2)  ~  2  •  (2.81 ) 

This  result  should  provide  a  reasonable  estimate  of  /,,,n_2)  to  within  a  factor 
of  sr  2.5.  The  expression  is  not  useful  in  consideration  of  specific  modes,  but 
will  provide  an  order  of  magnitude  view  of  7(,(„_2).  and  therefore  i l±(r=.2)i  for  a 
particular  range  of  I0  and  n  at  the  fundamental. 

An  example  of  this  scaling  is  graphically  displayed  in  Figure  2.5.  The  two 
graphs  in  the  figure  are  from  Gaponov’s  paper  [45]  describing  the  1975  Nusinovich 
and  Errn  results.  The  plot  at  the  top  is  applicable  to  the  fundamental  and  the 
bottom  graph  is  for  the  second  harmonic.  The  horizontal  scales  are  the  same 
between  the  two  diagrams,  but  the  vertical  axes  of  /,  have  different  log  scales. 

The  dashed  line  on  the  n  1  figure  follows  the  first  ridge  of  contours  rep¬ 
resenting  /„  optimized  with  respect  to  ?/  over  that  range  of  ft.  This  would  be  a 
typical  range  of  investigation  for  design  of  a  gvrotron  operating  at  the  fundamen¬ 
tal.  The  long  arrows  down  to  the  n  -  2  plot  show  how  specific  points  on  the  dashed 
curve  map  to  the  second  harmonic  by  the  derived  relations  (2.81)  and  (2.75).  For 
this  example,  (1S  =  0.38  was  used  in  the  calculation  of  /(j(„  2)-  The  dashed  line 


on  the  n  —  2  diagram  is  the  map  of  the  dashed  line  above.  The  vertical  bars 
extending  from  the  mapped  points  represent  the  estimated  range  of  error. 

An  interesting  aspect  of  the  example  presented  is  the  fairly  constant  value 
of  77  on  the  second  harmonic  dashed  line.  The  range  of  extended  from  45%  to 
72 %  at  the  fundamental,  yet  the  scaled  values  of  Jc,  and  [t  at  the  second  harmonic 
result  in  ?/ .  ~  38/7  t  2*7.  Even  including  the  rather  large  error  bars,  the  scaled 
values  still  remain  in  a  high  efficiency  zone  with  25%  <  rj_  •  50%.  This  would 
be  an  area  of  concern  for  megawatt,  n  -•  1  designed  gyrotrons.  Such  high  levels  of 
second  harmonic  efficiency  could  damage  components  or  degrade  efficiency  at  the 
fundamental. 


Chapter  3 

Description  of  the  Experiment 

The  primary  goal  of  the  experimental  investigation  was  to  determine  the  pres¬ 
ence  of  harmonic  emission  from  the  MIT  gvrotron.  Secondary  objectives  included 
frequency  measurement  for  mode  identification,  starting  current  determination 
and  output  power  level  measurement. 

The  MIT  gvrotron  is  described  in  Section  1.2  of  this  paper.  The  0,3,1  res¬ 
onator  cavity  introduced  in  the  same  section  and  discussed  throughout  Chapter  2 
was  used  exclusively  in  the  experimental  investigation.  The  whispering  gallery 
cavity  has  undergone  initial  testing  at  the  fundamental,  but  the  harmonic  content 
has  not  been  studied  to  date. 

This  chapter  will  explain  the  procedures  used  in  obtaining  the  experimental 
results.  The  first  section  details  the  magnetic  field  scan  to  detect  harmonic  emis¬ 
sion.  Section  3.2  discusses  the  frequency  measurement  technique  which  allowed 
mode  identification.  The  third  section  describes  the  starting  current  determina¬ 
tion  and  the  chapter  concludes  by  explaining  the  power  measurement  techniques 


3.1  —  The  Magnetic  Field  Scan  to  Detect  Harmonic  Emission 

The  primary  objective  of  the  experimental  investigation  was  fulfilled  by  a  scan 
over  the  normal  operating  magnetic  field  range  to  detect  the  presence  of  second 
or  higher  harmonic  modes.  The  fundamental  emission  was  mapped  simultane¬ 
ously  over  the  same  field  range  to  correlate  results  with  the  possibility  of  mode 
competition  between  the  fundamental  and  harmonic  modes. 

Figure  3.1  diagrams  the  apparatus  used  during  the  magnetic  field  scan.  The 
radiation  is  coupled  out  of  the  gyrotron  through  one  meter  of  2.2  cm  diameter 
oversized  cylindrical  waveguide.  The  fundamental  and  second  harmonic  receivers 
depicted  in  the  figure  were  placed  in  the  far  field  radiation  pattern  approximately 
16  inches  from  the  end  of  the  output  waveguide.  The  receivers  were  aligned  in  the 
same  horizontal  piane  as  the  waveguide  and  placed  symmetrically  at  18°  off  axis 
to  coincide  with  the  major  lobes  of  the  far  field  pattern. 

The  second  harmonic  receiver  consisted  of  a  25  dB  gain  VVR-3  band  horn,  a 
one  inch  section  of  WR-3  band  waveguide,  a  variable  dial  attenuator  and  a  video 
diode  for  a  detector.  The  operational  frequency  range  for  WR-3  band  is  220-325 
(ill/  with  a  cutoff  frequenct  of  172  CIlz.  The  diode  was  a  Baytron  model  4-06- 
005  x  wafer  in  a  Baytron  4-06- 105 A  detector  mount.  The  fundamental  receiver 
used  a  25  dB  gain  WR-0  band  horn,  a  dial  attenuator,  Hughes  wavemeter  and 
a  Hughes  model  47328H-11I1  video  diode.  The  WR-G  band  has  an  operational 
frequency  range  of  1 10  170  (Hlz. 

The  diode’s  were  chosen  to  be  optimized  for  (lie  frequency  band,  but  both 
showed  some  response  in  the  other’s  frequency  range.  Precautions  were  taken  to 
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ensure  only  second  harmonic  signals  entered  the  2uc  diode  and  the  fundamental 
diode  received  only  fundamental  emission.  In  the  first  case,  the  fundamental 
radiation  was  below  the  cutoff  frequency  for  the  WR-3  band  waveguide  and  dial 
attenuator  used  in  the  2 uic  receiver.  To  prevent  the  higher  harmonic  radiation  from 
entering  the  fundamental  receiver,  a  2.2  inch  thick  piece  of  plexiglass  absorber 
was  placed  in  front  of  the  receiver  horn.  The  attenuation  at  the  fundamental 
was  measured  by  comparing  relative  diode  signal  strengths  with  and  without  the 
plexiglass.  The  results  showed  an  8.0  dB  reduction  in  signal  at  128  GHz.  This  is  in 
agreement  with  the  i'~  scaling  of  the  absorption  coefficient  which  will  be  discussed 
in  Section  3.4.  The  success  at  preventing  leakage  of  the  higher  frequency  2 u)c 
radiation  was  indicated  by  placing  the  plexiglass  in  front  of  the  second  harmonic 
receiver.  This  resulted  in  the  complete  loss  of  the  strongest  second  harmonic  signal 
at  241  GHz,  representing  an  attenuation  of  over  27  dB.  As  a  final  precaution,  both 
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detectors  were  encased  in  a  layer  of  Eccosorb  material  to  prevent  reflected  RF 
signals  from  entering  the  diodes. 

The  dial  attenuator  for  the  second  harmonic  receiver  system  was  calibrated 
allowing  conversion  of  dial  settings  to  relative  signal  loss.  The  harmonic  emission 
from  a  Hughes  140  GHz  IMI’ATT  device  provided  the  radiation  during  the  cali¬ 
bration.  The  change  in  signal  st  rengt  h  as  a  fund  ion  of  dial  attenuator  setting  was 
recorded.  The  variable  attenuator  used  in  the  fundamental  receiver  system  w'as 
precalibrated  for  direct  dB  readings.  Throughout  the  magnetic  field  scan,  the  rel¬ 
ative  signal  strength  was  determined  by  adjusting  the  calibrated  dial  attenuators 
to  maintain  constant  diode  output  voltages.  This  method  was  preferred  since  it 
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precluded  damage  or  inaccurate  readings  caused  by  saturation  beyond  the  linear 
response  region  of  the  diodes. 

The  diode  signals  were  monitored  simultaneously  on  separate  channels  of  a 
dual  trace  oscilloscope.  Figure  3.2  shows  pulse  shapes  from  the  second  harmonic 
and  fundamental  diodes.  The  0. 8-1.0  n sec  pulse  lengths  indicated  in  the  figure 
were  typical  throughout  the  experimentation.  The  high  frequency  noise  prior  to 
the  pulse  is  pickup  from  the  high  voltage  power  supply.  The  repetition  rate  was 
*4.16  pulses  per  second.  The  cathode  voltage  Vc  =  64.3  kV  and  beam  current  7  =  5 
amps  listed  in  Figure  3.2  were  held  constant  through  the  range  of  magnetic  field 
values. 

The  magnetic  field  was  scanned  from  4.0  to  6.0  tesla.  The  lower  value  of  4.0 
tesla  was  chosen  to  ensure  adequate  beam  quality.  The  upper  limit  ensured  that 
the  fundamental  mode  frequencies  remained  below  cutoff  for  the  WR-3  band 
waveguide  used  in  the  2u;c  receiver.  The  magnetic  field  was  varied  in  0.044  tesla 
increments.  At  each  point  the  gun  coils  were  adjusted  to  optimize  the  signal  from 
the  second  harmonic  diode.  The  dial  attenuators  were  varied  to  keep  both  diode 
voltages  between  10-40  mV.  The  measured  signal  voltage  and  dial  attenuator 
reading  were  recorded. 

Figure  3.3a  graphs  the  results  from  the  fundamental  diode.  The  TE2,3.i  mode 
at  5.26  tesla  provided  the  zero  dB  reference  for  the  relative  diode  signal  plotted 
on  the  vertical  axis.  As  noted  on  the  figure,  the  diode  response  is  frequency 
dependent.  Therefore,  the  peaks  plotted  on  Figure  3.3a  do  not  provide  a  reliable 
measurement  of  mode  intensity.  For  example,  comparable  power  is  obtained  in 


Fig  3.2  —  Fundamental  and  Second  Harmonic  Diode  Pulse  Shapes 

1  —  5  amps,  Vc  =  64.3  kV,  B  =  5.26  tesla 


Horizontal  Time  Scale  —  0.5  psec  cm 


Top  Trace  -  2~v  Diode  Signal  into  280ft  Terminator 
\'crtical  Scale  20  mV /cm  with  3.3  dB  Attenuation 
TE4  51  Mode  at  263.31  GHz 


Bottom  Trace  —  wc  Diode  Signal  into  200ft  Terminator 
\ertical  Scale  -  10  mV  cm  with  44.6  dB  Attenuation 
7  F 2.3.1  Mode  at  137.35  GHz 


Fig  3.3b  —  Second  Harmonic  Emission 


Magnetic  Field  (Tesla) 


the  TEo.3,1  and  TE4)2,i  modes,  yet  the  latter  yields  a  significantly  lower  diode 
signal.  The  areas  on  the  graph  where  there  are  no  peaks  correspond  to  no  diode 
signal  with  the  dial  attenuation  reduced  to  zero. 

The  relative  diode  signal  from  the  second  harmonic  receiver  is  shown  in  Fig¬ 
ure  3.3b.  The  TE  11,2.1  mode  at  4.70  tesla  provided  the  zero  dB  reference  in  this 
case.  The  diode  response  is  frequency  dependent  so  conclusions  concerning  output 
power  should  not  be  drawn  by  the  relative  peak  heights.  This  is  particularly  true 
of  the  TE;-,,2.i  and  TE7.3,]  modes  at  4.12  and  4.40  tesla  respectively.  The  TE7  3J 
mode  at  226  GHz  is  on  the  edge  of  the  220-325  GHz  frequency  band  for  the  2wc 
diode.  The  TEc,,2,i  mode  at  209  GHz  is  actually  out  of  this  frequency  range  and 
the  diode  response  characteristics  could  be  quite  different  from  its  behavior  at  241 
GHz  for  the  TEn,2,i  mode.  The  output  power  measurements  using  a  calorimeter 
described  in  Section  3.4  verify  the  misleading  nature  of  the  diode  signal  by  show¬ 
ing  the  TEn.2.1  to  be  a  stronger  mode  than  the  TEg,2,i-  The  rest  of  the  peaks 
plotted  above  5.0  tesla  in  Figure  3.3b  are  within  the  operational  frequency  band 
of  the  diode  and  should  provide  some  measure  of  relative  mode  strength. 

The  modes  identified  on  both  Figures  3.3a  and  3.3b  were  determined  by  fre¬ 
quency  measurements,  for  the  more  stable  signals  al  the  fundamental,  the  Hughes 
wavemeter  included  in  the  fundamental  receiver  provided  an  adequate  frequency 
measurement.  However,  an  harmonic  mixing  system  yielded  higher  accuracy  and 
was  more  useful  for  some  of  the  weaker  or  more  erratic  signals.  This  frequency 
measurement  system  is  the  subject  of  the  next  section. 


3.2  —  Frequency  Determination  by  Harmonic  Mixing 

A  highly  accurate  and  sensitive  method  of  frequency  measurement  using  an 
harmonic  mixing  technique  was  employed  during  the  magnetic  field  scan  to  identify 
specific  second  harmonic  modes.  The  method  had  been  used  previously  for  real¬ 
time  spectral  analysis  of  far  infrared  laser  signals  |51j  and  has  also  been  used 
to  measure  the  bandwidth  and  frequency  pulling  of  the  MIT  gyrotron  at  the 
fundamental  52  . 

The  components  of  the  mixer  system  are  displayed  schematically  in  Fig¬ 
ure  3.4.  The  RF  from  the  gyrotron  is  heterodyned  in  the  harmonic  mixer  with  a 
harmonic  of  the  lower  frequency  local  oscillator  (LO).  The  resultant  intermediate 
frequency  (IF)  is  given  by  the  following  equations. 

t'ur  r—  upper)  ~~  k'lF  (3.1a) 

t'nr  —  ft  L' '^L1  ■  (lower)  ~  ^if  (3.1b) 

Here  . ( ,ippt-r i  and  */i,..(i..wer)  represent  the  fundamental  LO  frequencies  corre¬ 
sponding  1o  the  tipper  and  lower  sidebands.  The  IF  signal  is  amplified  and  then 
processed  by  a  surface  acoustic  wave  (SAW  )  dispersive  delay.  This  allows  the  IF 
frequency  to  be  displayed  as  a  function  of  time  on  an  oscilloscope. 

When  a  second  harmonic  frequency  measurement  was  desired,  the  2~y  diode 
signal  was  optimized  by  adjusting  the  magnetic  field  and  gun  coils.  The  receiver  for 
the  harmonic  mixing  system  was  positioned  in  place  of  the  fundamental  receiver 
shown  in  Figure  3.1  and  the  plexiglass  absorber  was  removed.  The  harmonic 
mixing  receiver  system  consisted  of  a  25  dB  gain  WR-3  band  horn,  a  two  inch 


section  of  WR-3  band  waveguide,  a  dial  attenuator  similar  to  the  :>ne  used  on  the 
2u;t.  receiver,  a  waveguide  transistion  from  WR-3  to  WR-6  and  a  Hughes  mode] 
47448H-1002  harmonic  mixer.  The  WR-3  waveguide  prevented  the  gvrotron's 
fundamental  emission  from  entering  the  harmonic  mixer.  The  effect  of  reflected 
RF  radiation  was  minimized  by  a  layer  of  Eccosorb  around  the  mixer. 

The  LO  signal  was  provided  by  a  YIG  luned  Gunn  oscillator  which  produced 
20  30  mW  of  output  power  in  the  12  16  GHz  frequency  band.  The  RF  frequency 
range  of  200  300  GHz  required  mixing  with  high  harmonics  of  the  local  oscillator, 
with  15  ?/,  ,  20.  The  fundamental  LO  frequency  was  monilored  by  an  E1P 

model  548  frequency  counter. 

The  variable  gate  coupled  through  an  RF  switch  on  Figure  3.4  allowed  an 
ad  justable  250  300  nsec  segment  of  the  gyrotron  pulse  to  be  isolated.  The  variable 

delay  permitted  various  portions  of  the  pulse  to  be  gated.  For  the  frequency 
measurements  during  the  present  investigation,  the  gate  was  set  to  a  300  nsec 
segment  coincidental  with  the  flattop  of  the  gyrotron  pulse. 

The  SAW  device  was  |,i\l>():  dispersive  delay  which  analyzed  IF  signals  in 
the  170  370  Mil/  range  by  delaying  them  as  a  linear  function  of  frequency.  The 
resultant  time  delays  ranged  Irom  1  to  20.3  p sec.  My  adjusting  the  horizontal 
time  scale  on  the  oscilloscope-  to  inc  lude  the  tolal  16.3  //sec  dispersive  delay,  the 
100  Mil/  SAW  bandwidth  was  display  ed.  In  practice,  a  2  p sec  div  time  scale  was 
converted  to  a  12.3  MHz  div  frequency  scale. 

Adjustments  in  the  LO  frequency  resulted  in  linear  changes  in  the  IF  fre¬ 
quency  by  equation  (3.1).  The  time  scale  on  the  scope  was  calibrated  by  changing 
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the  LO  frequency  until  the  displayed  IF  pulse  disappeared  symmetrically  about 
the  midpoint  of  the  horizontal  axis.  The  pulse  disappearance  was  due  to  exceeding 
the  470-370  MHz  bandwidth  limits  on  the  SAW  device,  thus  the  midpoint  was 
known  to  represent  an  IF  frequency  of  420  MHz.  Because  the  higher  frequencies 
received  less  delay  in  the  SAW  .  the  left  side  of  the  scope  represented  uxv  420 
MHz  and  the  right  side  of  the  scope  displayed  signals  with  120  MHz. 

The  procedure  for  frequence  determination  involved  scanning  the  LO  fre¬ 
quency  until  an  IF  pulse  appeared  on  tin  -cope.  The  puKe  \\ a-  centered  and  the 
LO  frequency  was  recorded.  The  pulse  was  identified  a-  an  upper  or  lower  side¬ 
band  by  increasing  the  LO  frequence  slight le.  If  the  pulse  moved  right  indicating 
a  decrease  in  IF  frequency,  it  was  identified  as  tin'  lower  sideband  and  the  plus 
sign  in  equation  (3.1b)  applied.  If  the  pulse  moved  left  when  the  LO  frequency 
was  increased,  it  was  the  upper  sideband  and  the  minus  sign  in  equation  (3.1a) 
was  appropriate. 

The  local  oscillator  haiinonu  nutubet  | n.  )  could  easilx  be  determined  be- 
«  au-c  i'  uas  ah\a\-  adi'i-ted  s,i<  |,  :|1(l;  r  120  MHz.  By  locating  both  upper 
and  lower  sideband  .  "  <  nuh’  In  (  .m  ula'<  d  !>\ 
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(3.2) 


At  this  point,  even  thing  on  tin  right  side  of  equations  (3.1a)  and  (3.1b)  was 
determined  allowing  calculation  of  the  111  frequency  In  either  expression.  The 
instrumental  limit  on  the  arcurace  of  t  lie  frequence  measurement  was  approxi- 
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of  the  lack  of  significant  simultaneous  fundamental  oscillation.  For  the  TEj  1,2,1. 
the  second  harmonic  power  accounted  for  90-95%  of  the  total  output  power.  This 
percentage  was  slightly  lower  in  the  case  of  the  TE11.2.1.  with  78%  attributable  to 
the  second  harmonic  mode. 

Considering  the  range  of  values  listed  in  Table  3.2.  the  output  power  claim  of 
k\\  at  211  CM/  in  the  TE 1 1/...1  is  justified.  It  is  this  value  which  appears  on 
f  igure  I  lh  showing  gvrotron  performance  al  harmonics.  This  result  at  6  amps 
« i»rre*|M»t»d*  to  .»n  overall  efficiency  r/  s:  (i.T’V.  It  i"  unknown  if  the  slightly  lower 
output  |«*«wrr»  Mu  a*ufe«|  at  7  \  for  the  TEu.2.1  are  due  to  day-to-day  variations 
in  gvrotron  performance  or  a  phv-iral  effort  such  as  mode  suppression  at  higher 


rurrrnt  value*  1  hr  T I,  .  1  output  |Htwer  value  of  *s  15  kVV  translates  to  an 
overall  rfliornrv  of  approxirmilelv  l.Vi . 


lower  than  expected  efficiencies. 

Because  the  theoretical  predictions  of  Chapter  2  considered  only  the  resonalor 


geometry,  the  problem  of  radiation  loss  or  absorption  at  the  output  window  or 


in  coupling  with  the  output  waveguide  was  not  considered.  The  quartz  output 
window  was  specifically  designed  for  the  140  GHz  radiation  of  the  fundamental.  It 
is  unknown  exactly  how  much  increased  reflection  or  absorption  may  be  occurring 
in  the  window  at  the  higher  second  harmonic  frequencies.  Increased  reflection 
would  raise  the  Q  and  change  the  theoretical  calculations  and  increased  absorption 
would  obviously  lower  the  measured  radiation  levels. 

The  guiding  center  radius  r.  and  ratio  3  are  parameters  which  enter  into 
both  the  linear  and  nonlinear  theory  calculations,  yet  they  were  not  measured 
directly  in  this  experiment.  The  values  which  were  used  in  the  starting  current 
and  efficiency  calculations  were  design  figures  based  on  electron  gun  computer 
simulations  and  adiabatic  theory.  As  an  example  of  how  a  reasonable  change 
in  these  parameters  could  affect  the  differences  between  theory  and  experiment, 
equation  (2.32)  shows  that  Lr  ~  (/?_/;?  )"a  for  the  second  harmonic.  With 
ihe  design  value  of  . 1  d  -  1.49.  the  TEi  1.2.1  mode  has  a  predicted  value  of 
:  0.00  A.  However,  if  the  actual  ratio  of  ;  3  -  2.0.  the  change  would  yield 

/  T  0.37  A  which  is  in  much  better  agreement  with  the  experimentally  measured 
starting  current  of  0.2  A  for  the  TEj  1.2.1  inode.  A  consideration  of  this  sort  will 
be  analyzed  in  Section  4.1  concerning  the  value  of  the  guiding  center  radius  rf. 

This  example  is  a  simple  illustration  of  how  a  change  in  one  parameter  value 
could  account  for  some  of  the  difference  between  theory  and  experiment  in  the 
starting  current  for  one  isolated  mode.  In  general,  the  increase  in  3 .  /(tv  men¬ 
tioned  would  prove  counter  to  the  overall  observation  of  lower  than  predicted 
levels  of  second  harmonic  emission  and  reduced  efficiencies.  This  chapter  will 


analyze  physical  reasons  which  could  explain  the  general  trend.  The  discussion 
will  focus  on  two  separate  phenomena  as  possible  explanations.  The  first  will  be 
beam  radius  effects  including  the  actual  beam  placement  and  the  consideration  of 
the  finite  beam  width  in  calculations.  The  second  topic,  which  will  be  discussed 
in  Section  4.2.  is  the  possibility  of  mode  competition  between  fundamental  and 
second  harmonic  modes  as  well  as  between  neighboring  second  harmonic  modes. 

4.1  Finite  Beam  Width  Effects 

Kxaminalion  of  Table  3.1  listing  the  observed  second  harmonic  modes  reveals 
a  trend  in  the  radial  mode  index  p.  With  the  exception  of  the  TE^.s.i.  each 
mode  identified  had  radial  mode  numbers  of  p  2  or  p  3.  The  linear  theory 
analysis  forecast  modes  above  threshold  wilh  2  p  7.  and  predicted  the  modes 
with  p  -  4  to  have  the  lowest  starting  currents.  This  discrepancy  prompted 
consideration  that  the  beam  radius  was  possibly  different  from  the  value  used  in 
the  theoretical  work,  or  some  other  interaction  was  occurring  based  on  the  radial 
thickness  of  the  beam  which  would  suppress  the  higher  order  radial  modes. 

The  guiding  center  radius  (r, )  enters  into  both  the  linear  thcon  starting 
current  calculations  from  equation  (2.321  and  the  nonlinear  tlicon  calculations 
of  /  in  equation  (2.(»8)  by  the  same  G  factor  defined  in  equation  (2.28).  This 
term  represents  the  efficiency  of  interaction  between  the  electron  beam  and  the 
I? I"  fields.  It  is  actually  only  the  numerator  of  (7  whirl)  measures  the  strength 
of  interaction  by  G  —  (A-  r,).  with  the  denominator  representing  iho  mode 

stored  energy.  Since  /. T  ~  G  1  and  /,  -  G.  an  increase  in  G  lowers  the  starting 
currents  and  increases  /,..  This  increase  in  results  in  an  increase  in  t  he  predicted 

KM) 


*•  .  »  -  »  .  * 


value  of  r/  throughout  the  range  of  physically  realizable  beam  current*  in  the 
MIT  gyrotron.  The  discussion  in  this  section  will  consider  the  physics  which 
could  account  for  an  effective  lowering  of  the  C  factor,  thereby  raising  the  starting 
currents  and  reducing  the  predicted  second  harmonic  emission.  Although  the  focu* 
will  be  on  raising  the  starting  currents,  the  same  arguments  will  result  in  lower 
theoretical  efficiencies  by  the  /  parameter  lined  in  the  slow  time  scale  theory. 

The  minimum  beam  radial  thickness  is  two  I, armor  radii  (2r,  ).  Tin1  actual 
t  hickness  is  greater  due  to  the  finite  size  of  the-  electron  emitter  sttip.  space  c  barge 
and  other  effects.  The  simplification  which  allowed  determination  of  by  analy  tic 
methods  was  that  the  interaction  could  he  represented  using  an  infinitely  thin 
annular  beam  with  radius  r,.  so  it  is  reasonable  to  assume  that  some  physics 
considerations  were  lost  by  discarding  the  radial  beam  thickness.  For  the  MIT 
gyrotron.  the  electron  beam  width  (thickness}  predicted  by  computer  calculations 
is  3.o  r. .  Since  r,  r  .  it  is  evident  that  the  beam  thickness  was  changing 
throughout  the  experimental  magnetic  field  scan.  I’sing  tin-  field  values  li*ied  in 
Table  2.2a.  the  beam  widths  range  from  0.003  cm  for  tin  TK|  •  t  mode  at  1.05-1  T 
t“  0.012  cm  lor  the  Ti.|*..j,i  mode  at  0.017  T.  An  average  value  (or  the  beam 
width  ot  approximately  0.05  cm  represent*  ^  25' 7  ol  tin-  r.  0.107  cm  value-  used 
in  tin’  theoretical  calculations. 

Figure  1.1  is  representative  of  the  extent  of  variation  in  calculated  starting 
current  which  would  result  by  using  different  value’s  of  r,  within  the  range  of  beam 
thickness.  The  two  modes  represented  are  the  two  senses  of  op|>osite  azimuthal 
rotation  for  llie  TK.|..|.|.  The  presented  values  of  /  T  are  valid  at  li  1.350  T  as 


|r»trd  in  Tab le  2  2*  The  brain  hmIiN  at  thr*  magnet*  hrld  ulor  »►  O.OW  cm.  *u 
l he  edge*  ol  the  plotted  line*  at  U  lb»  cm  *uii  0  J2I*  cm  represent  tbe  mint  and 
oulrr  icdiiilrtiir*  ol  (hr  br«m  Thr  ilcpli  turiilrtrij  uii  (hr  t*lur  d  f ,  U  IU?  t  tu 
mh*h  **a*  u*ct)  in  Ihr  Marling  currmi  and  Htinruii  ««kuU(w)t»  in  (  hnjnn  2 
The  2  1  imp  value  of  J.,  lot  the  TI-*  «  i  at  ••  H*7  «m  i»  hgui*  Jotrd  in 
T«l>ir  2  2*  The  da*hed  hnnionul  line  irjiirwun  ll«e  huitmu  ’•  amp  oper  at  mjj 
current  during  the  experimental  imr»iig*iii«i 

It  appear*  ol<t  nm<  lrt*m  figure  I  I  that  vet*  huh  >idoiiiiat*cu  uiim  rtiung 
the  Iteam  ititetactNiiii  i<  g  waned  from  a  *impk  «ahul*i**n  of  /  (•*>*«(  «ti  ih» 

ndinitelv  thin  line  at  ll  |"*7  cm  The  plot*  extending  above  th«  graphed  limn  «d 
/  KNNI  atlijr*  air  acluallv  inltlOt*"*  i«««ri|  In  the  \#lgr»  ol  f,  «ln(h  tc'Mill 

m  J;,  ,,|*  r.)  0  m  the  denominator  \  minor  increase  in  t,  d  »  -0  011“  «m 

mould  place  the  calculated  value  o(  /  ,  above  the  '*  amp  line  for  the  Th«  t  k,  >  The 
previous  de*f  u**'wn  concerning  the  degree  erf  Know  ledge  about  the  actual  value  d  t , 
i»  now  relevant  *ince  minor  change*  re*ult  in  «u»h  major  ddfetrucr*  jn  cahuiaicd 
value*  of  /  . 

The  diagram  *uggeMrd  two  po*«iblc  improv eluent*  t..  Mu  value  ol  <«  in  Ik 
inr  or|K»r«led  into  the  Matting  »  uflrnt  and  rfliciencv  calculation*  Tin  ht»t  *tcp 
wa*  *implv  an  attempt  at  improving  the  acrurarv  ol  the  a*«uuwd  value  of  r  l.aod 
on  the  actual  experimental  magnetic  field  value*.  The  «eeond  approach  involved 
taking  the  integral  average  of  the  CV  factor  over  the  finite  beam  vvidth  to  account 
for  the  change*  in  the  interaction  strength  with  different  segment*  of  the  Imam 
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This  limited  apparent  success  is  marred  by  the  overall  result  that  the  finite 
beam  thickness  factor  did  not  eliminate  any  modes  by  increasing  their  theoretical 
starting  currents  over  the  5  amp  experimental  operating  current.  It  actually  added 
the  TKi.t  .n  i  ntode  to  the  list  of  those  predicted  to  be  above  threshold.  In  simple 
terms,  the  beam  width  averaging  increased  the  starting  current  for  the  low  starting 
current  modes,  but  decreased  /  ,  for  the  high  starting  current  modes.  The  range 
of  values  of  0.  II  A  /  ;  'ifXX  1  A  for  the  infinitely  thin  beam  calculations  was 
compressed  to  0  I*  \  /  310  A  with  the  beam  width  factor.  Additionally, 

the  factor  did  not  afvcavs  preferentialic  lower  starting  currents  of  the  low  radial 
modes  and  raise  /  ..  for  tin  higher  order  modes,  as  evidenced  bv  tliefXi'V  reduction 
in  I  «  for  the  Tf.t  ,  Thi*  mode  will  serve  as  a  useful  summary  for  the  beam 
radio*  considerations  of  this  section  and  provide  some  insight  into  the  apparent 

fadings 

figure  I  :s  depict*  the  predicted  starting  rurrent  for  the  TK » .i  modes  at 
1  *r-l  T  The  left  and  tight  edges  of  the  curves  represent  the  inner  ami  outer  radial 
edge*  of  the  3  *.  *  lieamwidlh  at  t|o«  magnetic  field  calm  The  initial  calculation 
of  l  ..  for  thr«  roo«ir  in  Chapter  2  produce*!  /  *•  <•  \  {or  ilu  I  l.i,  ]  .  Thi* 

vain*  vea*  t»a«cd  on  a  '“thin'*  l>«  am  mil  * .  O.I'lT  « m  ami  i»  ide  milled  on  I  rgiife  •!..'{ 
w  >rl  'hr  -  I  Th*  ifaitinj  oinrisi  of  < f»»  Tl.|  i  a*  thi»  calm  of  t,  i*  tr  1000  A. 
indicated  b\  rh*  2  on  the  grapn  The  improved  rsiimatc  of  r  O.lfll  cm  still 
shnwrd  the  Tl.  (  a'  the  Sower  starting  current  mode,  hut  the  value  had  more 

than  doubted  lo  /  |.t  \  which  i‘  depictrd  fiy  the  ?■  .1  The  reduction  in  r 

from  0  107  err*  to  ft.  |*l|  no  resulted  »*■  /  -  for  the  Tf,  ,  rnoele  increasing  from 
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1000  A  —  oc,  then  decreasing  back  down  to  the  value  of  7.T  56  A,  shown  by 

the  #4  in  the  diagram.  Finally,  the  beam  width  averaged  starting  current  values 
showed  the  two  modes  approximately  equal  with  7.T  ~  2.3  A  for  the  TK j .... i ( -  , 
and  7.T  -  2.2  A  for  the  TEi>(>,k_).  Although  the  central  region  of  the  beam 
interacts  very  weakly  with  the  fields  as  evidenced  by  the  high  “thin”  beam  7.T 
values,  the  inner  and  outer  edges  of  the  !  cam  are  positioned  to  interact  favorably, 
thus  bringing  the  lull  beam  width  averaged  starting  current  down  to  2.2  A. 

Since  the  two  curves  for  the  rotating  modes  art'  almost  symmetric  about  the 
central  value  of  r,  0.191  cm.  it  is  reasonable  to  assume  that  their  interaction 
with  the  fields  over  the  full  beam  width  would  be  similar.  The  beam  width  aver¬ 
aged  starting  current  calculations  support  this,  whereas  the  delta  function  view  of 
the  beam  radius  yields  orders  of  magnitude  differences  in  the  predicted  starting 
currents  for  the  twfo  senses  of  rotation.  Jf  is  felt  that  the  averaged  values  are  more 
physically  realistic,  even  though  the  TKi.r.j  mode  was  not  observed  experimen¬ 
tal^.  The  lack  of  mode  excitation  is  probably  due  to  phenomena  unrelated  to 
beam  radius  effects. 

A  possibility  not  considered  in  t  he  t  lieoret  ical  calculal  ions  is  i  hat  t  here  may  be 
standing  wave  st  rucl  ures  in  t  he  cavity  as  opposed  to  t  Ik  total  mg  waves  assumed  for 
the  nnna/.imuthalh  symmetric  modes.  Thi-  would  lower  the  predicted  interaction 
because  of  electromagnetic  field  nulls  in  the  azimuthal  direction  which  do  not 
exist  if  the  wave  is  rotating.  Experimental  measurements  of  the  far  field  pattern 
indicate  that  there  are  rotating  modes,  but  ibis  does  not  preclude  the  possibility 
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of  standing  wave  modes. 


Although  possible  improvements  to  the  averaging  technique,  such  as  including 
a  spatial  distribution  function  for  the  electrons,  might  yield  minor  improvements, 
the  limited  success  at  predicting  experimental  results  suggests  consideration  of  dif¬ 
ferent  explanations.  The  degree  of  success  of  the  finite  beam  width  considerations 
is  limited  to  the  range  of  applicability  for  the  theories  covered  in  Chapter  2.  These 
were  single  mode  analyses,  so  the1  failure  of  the  beam  width  corrections  to  account 
for  the  experimental  observations  is  possibly  due  to  the  failure  of  the  theory  itself, 
if  multimoding  or  mode  competition  effects  were  present.  These  possibilities  are 
discussed  in  the  following  section. 

4.2  Mode  Competition 

Theoretical  analyses  of  multimoding  and  mode  competition  in  gyrotrons  have 
been  conducted  previously,  such  as  the  investigations  by  C.  S.  Nusinovich  55  and 
Kreischer  et  a!.  5(>  .  These  studies  have  shown  thai  the  phenomena  can  have  sig¬ 
nificant  impact  on  the  mode  oscillation  characteristics  of  gyrotrons.  Nusinovich 
cites  three  primary  manifestations  ol  multimoding.  these  being  mode  suppression, 
mode  enhaiu  ement  and  mod*  locking.  Ol  t lit*  three,  mode  suppression  is  the 
dominant  i  lb  <  :  nu  mmii  •  ha\t  <  io-i  o-)  illation  Itequeiu  ies.  This  is  char¬ 

acterized  by  an  iKii\«i\  oscillating  mode  limit  ing  t  lie  oscillal  ion  of  nearby  modes 
which  might  oth<  tvCm  be  above  threshold.  The  predictions  of  these  events  are 
clearly  outside  tin  realm  of  the  theories  presented  in  Chapter  2.  In  a  similar  man¬ 
ner.  1  lie  existence  of  mode  competition  effects  could  severely  reduce  the  accuracy 
of  the  calc ulat ions  based  on  those  theories. 


The  study  by  Kreischer  and  coworkers  cited  above,  included  an  experimental 
investigation  of  multimoding  at  the  fundamental  with  the  same  0.3.1  cavity  used 
in  this  study.  There  were  large  areas  of  mode  overlap  which  were  predicted  by 
theory  and  supported  by  the  experimental  results.  Additionally,  the  fact  that  the 
TE,..8.i  mode  never  fully  acheived  the  high  efficiencies  otherwise  predicted  was 
attributable  to  its  oscillation  zone  overlap  with  the  TK;.:i,i  mode. 

In  the  present  study,  the  diode  signal  traces  such  as  those  shown  in  Figure  3.2. 
confirm  the  existence  of  multimode  oscillation  between  fundamental  and  second 
harmonic  modes.  Additionally,  the  simultaneous  frequency  measurement  of  the 
TK. i i  and  TKir, >2  ,  inodes  at  the  same  magnetic  field  value,  is  a  clear  example 
of  multimoding  between  second  harmonic  modes.  The  existence  of  this  multimode 
oscillation  points  to  the  possibility  of  mode  suppression  or  enhancement  between 
fundamental  and  second  harmonic  modes,  or  between  separate  second  harmonic 
modes. 

The  mode  competition  theory  developed  by  Nusinovich  is  valid  for  arbitrary 
harmonics,  but  analytic  methods  for  applying  the  theory  have  only  been  developed 
at  the  fundamental  with  some  simplifying  assumptions,  such  as  a  tiniform  distri¬ 
bution  for  the  axial  field  profile.  Quantative  methods  are  therefore  not  available 
for  a  detailed  analysis  of  the  effects  of  multimoding  at  the  second  harmonic  or 
mode  competition  between  fundamental  and  second  harmonic  modes.  However, 
qualitative  arguments  will  be  presented  as  well  as  the  experimental  evidence. 

In  a  magnetic  field  region  which  supports  both  fundamental  and  second  har¬ 
monic  modes,  the  fundamental  will  generally  begin  oscillation  prior  to  the  second 


harmonic  because  it  has  higher  gain,  a  lower  starling  current  and  a  wider  excitation 
range.  This  has  been  experimentally  confirmed  An  example  of  this  phenomenon 
is  shown  in  Figure 3.2.  where  the  fundamental  TK;  j  »  diode  signal  turns  on  earlier 
and  has  a  longer  pulse  length  than  the  second  harmonic  TFI4.ii  mode.  Although 
the  second  harmonic  mode  may  still  oscillate,  the  established  fundamental  mode 
will  tend  to  deplete  the  available  electrons  which  might  otherwise  exchange  en 
ergy  with  the  second  harmonic.  The  result  would  be  reduced  efficiency  and  output 
power  at  the  second  harmonic  from  that  which  would  be  predicted  in  tin  absence 
of  the  fundamental  mode. 

The  experimental  results  of  thi*  study  appear  to  siip|»orl  the-  idea,  in  pat- 
titular  the  case  of  the  fundamental  suppressing  the  second  harmonic  emission 
Figure  1. 1  reproduces  the  magnetic  field  scan  diode  signal  plots  from  Figures  3.3a 
and  3.3b.  Both  «•  and  2-\  diode  signals  are  shown  on  the  same  graph,  so  the 
vertical  scale  has  been  eliminated  to  avoid  misleading  estimates  of  relative  mode 
strength.  As  menlioned  in  Chapter  3.  relatively  few  com  fusions  may  In*  drawn 
concerning  output  power  by  the  diode  signal  even  from  the  same  diode,  vet  the 
TF,,  ..  ,.  TK  |  and  TKr  >  t  modes  provided  the  most  consistently  stable  diode- 
signals  at  the  second  harmonic.  The  higher  magnetic  field  2»,  modes.  TK  \ 
through  the  TF.,-. were  generally  characlerimf  by  erratic,  low  level  diode  sig¬ 
nals.  present  over  a  much  more  limited  range  of  magnet ir  field  and  gun  coil  settings. 
As  seen  in  Figure  4.4.  these  high  magnetic  field,  low  level  2~v  modes  occur  at  the 
same  field  values  as  significant  levels  of  fundamental  radiation,  in  particular  the 
Tlv  :  , .  TK  .  i  and  TK-.  ■<  \  modes  By  contrast,  the  consistently  stable  TK,, .2  .1. 
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magnetic  IkW  range.  flight  It  uurr  than  onr  half  af  all  poMiblr  third  btrawtitt 
mode*  would  hr  cm  dialing  *>  6*e  a  tup*  and  M3  IV  brant  idiatr 


The  fourth  harmonic  rnuin  ImH  to  Table  ’■»  2  *how  no  mode*  above  thrrdndd 
at  five  amp*  brant  current  Thr  minimum  »tatimg  currrot*  to  tht»  (#**•  are  on 
the  order  of  il  A  The  lari  of  ant  predated  cm  dial  ton  at  the  fourth  iamotin 
eliminated  the  need  lor  convaderalwtn  of  rttode*  «t«th  n  4 

Mtiornew*  were  calculated  lor  the  third  harmccn*  mode*  predarted  to  be 
alcove  threshold  utilig  the  *lra»  litter  a«lr  ncaNlinrat  I Wm t  njeatibcd  tn  t'haplef  2 
**I>cC(l»f alb,  the  *hoft  integration  length  c d  p  -/#»  wa*  u««rd  The* 

t  t 

nonlinear  thcon  generally  predated  tbr  hghrcl  rHaeMa*  w ben  compared  with 
the  long  integration  length*  and  the  fa*«  tm*  wtalr-  ihrata*  a  bo  u*ed  in  Chapter  2 
The  result*  of  the  third  barman*  rlhcirocv  cahulat am*  are  presented  in  Table  U 
lor  \  M  3  hV  and  /  %  \ 

Thr  third  harmonic  perpendicular  Hlicimcie*  (»»  I  were  *u»pmi«gb  high  lot 
the  mode*  con*»drrr«l  V  *ern  m  the  table.  21**  n  II**  lliwevn.  the  low 
ohmir  Q  factor*  relator  to  thr  diffractive  {/'*  fee*  the  n  3  mode*  *tgnibc  antic 
rriliicec)  the  nutpiil  rfheieOfur*  to  the  Iff  .  *,  V  .  tang*  Ofncotj*H  pc*  e*fh- 

rtrpoc*  could  If  ralrulatrd  for  fourth  ha'oxu.n  n«dc»  a'  (m  am|»*  «ipce  then 
wrrr  all  lie  low  th  re*  hold.  bowevrt  the  efbeerCHc  wa«  calculated  leer  the  TK,«  ‘  ti  - 
af  7  amp*  u*ing  the  «arrie  *|r»w  lime  *ealr  then**  Thr  predicted  value  of  p  «»• 
Ml'  I  re* tilling  in  an  overall  rBkrwt  n  3.3* * 
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The  lark  of  any  observed  third  harmonic  emission  is  probably  due  to  mode 
suppression  considerations  as  discussed  in  the  previous  chapter.  Both  the  fun¬ 
damental  and  second  harmonic  modes  have  higher  gain  and  should  suppress  the 
third  harmonic  oscillation.  The  modes  listed  in  the  Tables  5. 1-5. 3  should  occur 
at  ss  4.#  T.  but  this  coincides  with  the  second  harmonic  TE|  |,2.|  mode.  As  seen  in 
Figure  1.1.  there  are  no  magnetic  field  regions  between  1.0  6.0  tesla  where  there 
is  neither  fundamental  nor  second  harmonic  emission,  thereby  making  third  har¬ 
monic  (filiation  unlikely.  The  results  in  Tables  5.1  based  on  the  single  mode 
analysis  did  not  consider  these  mode  competition  possibilities. 


Th<*  present  study  has  confirmed  the  existence  of  second  harmonic  emission 
from  a  high  power,  high  frequency  gyrolron  designed  for  the  fundamental.  The 
second  harmonic  radiation  was  within  the  high  efficiency  zone,  resulting  in  approx¬ 
imately  2*>  k\\  of  output  power  at  24 J  GHz.  The  linear  and  nonlinear  theories 
which  were  applied  indicate  that  efficient  second  harmonic  emission  will  be  the 
general  result  rather  than  the  exception  in  overtnoded.  high  power  gyrotrons.  Fu¬ 
ture  megawatt  power  lex  el  gyrotrons.  designed  for  operation  ai  the  fundamental, 
will  have  to  consider  t he  effects  of  this  harmonh  oscillation  or  develop  methods 
for  suppressing  these  modes. 

The  harmonic  emission  observed  in  t  his  invest  igat  ion  was  identifiable  as  cavity 
modes.  That  is.  each  mode  frequency  _  =:  2~  could  be  predicted  based  on 
the  scaling  of  the  Hesse!  function  zero.  i',,,,..  The  distinction  here  is  between 
these  cavity  modes  and  parametric  effects  which  have  recently  been  identified 
in  the  radiation  from  a  2*  GHz.  cw  Yarian  gyrotron  to  be  used  for  ECRH  on 
the  TARA  tandem  mirror  experiment  at  MIT.  These  latter,  parametric  effects 


are  characterized  by  a  second  harmonic  frequency  measured  as  exactly  twice  the 
frequency  of  an  oscillating  fundamental  mode.  Experimental  work  is  continuing  to 
investigate  this  phenomenon,  and  separate  physics  issues  are  involved  than  with 
the  cavity  modes  found  in  the  gyrotron  used  for  this  study. 


6.1  —  Summary  and  Further  Conclusions 


The  theoretical  investigation  separately  considered  two  different  resonator 
cavities.  The  0.3.1  cavity  was  designed  for  scaling  up  to  megawatt  power  levels  and 
cw  applications.  The  relatively  small  beam  radius  to  wall  radius  ratio  naturally 
yielded  a  dense  mode  structure,  with  the  possibility  of  36  second  harmonic  modes 
between  4.0  and  6.0  tesla.  The  whispering  gallery  cavity  was  designed  for  reduced 
mode  competition  at  the  fundamental,  and  the  larger  beam  to  wall  radius  ratio 
allowed  a  relatively  sparse  spectrum,  with  only  half  the  possible  2uc  modes  as  the 
0.3.1  cavity  over  the  same  magnetic  field  range.  For  both  cavities,  the  CAVRF 
computer  code  predicted  the  4.0  6.0  tesla  second  harmonic  mode  frequencies  to 
be  in  the  203  303  CJIIz  range. 

A  linear  analy  sis  determined  t  he  t  hreshold  of  om  illat  ion  <  ondil ions  lor  t  he  2«*.\ 
mode'  in  both  cavities.  The  second  harnionit  spectrum  ol  ihe  0.3. 1  cavity  was 
predicted  to  be  heavily  overmoded.  with  28  of  the  36  possible  inodes  in  the  4.0-6.0 
tesla  range  having  starting  currents  of  less  than  five  amps.  The  whispering  gallery 
cavity  calculations  showed  a  greatly  reduced  mode  spectrum,  with  only  seven 
modes  predicted  to  be  above  threshold  at  a  beam  current  of  five  amps.  Similarly, 
the  minimum  theoretical  starting  currents  were  lower  for  the  0.3.1  cavity,  with 
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7<T  s»  0.5  A.  as  opposed  to  the  whispering  gallery  resonator  which  had  a  minimum 
1<T  ~  3.5  A. 

The  nonlinear  theory  was  employed  for  efficiency  calculations  for  all  the  sec¬ 
ond  harmonic  modes  predicted  to  be  above  threshold  by  the  linear  analysis.  A 
slow  time  scale  formulation  allowed  efficiency  calculations  from  computer  gener¬ 
ated  graphs,  based  on  the  assumption  of  a  Gaussian  longitudinal  field  profile. 
Application  of  the  slow  time  scale  theory  only  required  calculating  the  normalized 
operat  ing  current  /,,  and  normalized  cavity  lengt  h  //.  The  second  harmonic  modes 
predicted  to  be  above  threshold  in  the  0.3.1  cavity  had  total  efficiencies  ranging 
from  10  25*7.  The  whispering  gallery  cavity  slow  time  scale  theoretical  results 
were  in  the  5  10‘7  efficiency  zone.  A  fast  time  scale  computer  code  using  the 
cold  cavity  CAVRF  field  profile  was  also  employed  for  efficiency  calculations.  The 
resulting  efficiencies  were  lower  than  the  slow  time  scale  values  by  2  10  percentage 
points  for  both  the  0.3.1  and  whispering  gallery  cavities. 

Starting  currents  and  efficiencies  wen*  calculated  for  a  limited  number  of  third 
and  fourth  harmonic  modes  using  the  0.3.1  cavity.  The  lowest  third  harmonic 
starling  currents  were  approximately  2  A  and  the  fourth  harmonic  had  the  lowest 
values  of  7.  -  si  5.-1  A.  The  results  indicated  t  hat  slighl  ly  more  than  half  of  the  third 
harmonic  modes  between  -1.0  6.0  tesla  should  be  above  threshold  al  five  amps.  A 
very  small  percentage  of  the  fourth  harmonic  modes  should  be  excited  at  that 
beam  current.  Output  efficiencies  for  the  third  harmonic  modes  were  calculated 
to  be  in  the  5- 10*7  range  using  the  slow  time  scale  nonlinear  theory. 


The  final  theoretical  analysis  was  a  scaling  estimation  for  second  harmonic 
efficiency  based  on  operating  parameters  at  the  fundamental.  The  slow  time  scale 
efficiency  plots  provided  an  easy  method  of  mapping  regions  of  interest  at  the  fun¬ 
damental  to  the  appropriate  efficiency  zone  at  the  second  harmonic.  The  method 
supported  the  quantitative  results  of  the  0.3.1  cavity  by  predicting  high  efficien¬ 
cies  at  the  second  harmonic.  An  example  presented  for  a  device  optimized  at 
the  fundamental  showed  a  remarkably  constant,  high  level  of  second  harmonic 
efficiency,  even  to  the  point  where  the  fundamental  efficiency  had  decreased  by 
approximately  30'  <  .  'Phis  motivates  the  concern  over  harmonic  emission  for  low 
ohmic  loss,  low  Q,  overmoded  resonators  designed  for  fundamental  operation. 

The  experimental  investigation  was  designed  to  confirm  the  predicted  exis¬ 
tence  of  the  harmonic  emission  and  provide  some  test  of  the  theoretical  calculations 
for  the  0.3. 1  cavity,  A  magnetic  field  scan  over  the  4.0  6.0  tesla  range  at  64.3  kV 
and  with  '»  amps  of  beam  current  was  conducted,  using  video  diodes  to  detect  both 
fundamental  and  seeond  harmonic  emission.  Kight  2~  modes  w<  te  identified  by 
an  harrtionh  mixing  frequency  measurement  sy  strut.  with  the  frequencies  gener¬ 
ally  within  (».r.  of  the  podnied  *  AX’  I  >  I  value-.  This  mixing  system  proved  to 
be  invaluable  due  to  its  ;n  (  ur.u  \  ami  ,i i *i I i : v  i<>  measure  low-level,  erratic  signals. 

I  h«  starting  current  measured  ior  the  Tl,|t  ,  mode  was  0.2  A.  which  was  a 
factor  of  three  lower  than  the  theoretical  value .  and  of  the  same  order  as  /.T  for 
fundamental  modes.  Power  measurements  of  the  TK11.2.1  and  Tlv.  •>  \  modes  were 


performed  using  calorimeters.  The  reflectivity  of  the  calorimeters  at  these  frequen¬ 
cies  was  analyzed  using  dispersive  Fourier  transform  spectroscopy,  and  power  dis¬ 
crimination  between  the  fundamental  and  second  harmonic  was  accomplished  by 
using  various  thicknesses  of  plexiglass  absorber.  The  results  revealed  the  TEu.s.i 
mode  producing  ^  25  kW  of  output  power  at  241  GHz.  with  the  TKy.z.i  emitting 
15  k\N  at  209  GHz.  These  results  represent  the  highest  output  powers  at  high 
frequencx  (  100  GHz)  detected  in  a  gyrotron  operating  at  harmonics  above  the 

fundamental. 

The  detection  of  second  harmonic  emission  and  the  power  levels  achieved 
represent  the  most  significant  aspects  of  this  stud).  However,  less  than  one  third 
of  the  %  2~\  modes  predicted  to  be  excited  were  observeel  and  the  measured 
efficiencies  were  a  factor  of  three  lower  than  the  theoretical  values.  Additionally, 
no  higher  harmonic  (n  ■  2)  modes  were  observed  although  a  number  of  third  har¬ 
monic  modes  were  predicted  to  be  above  threshold.  Two  possible  explanations  for 
this  difference  were  investigated.  The  first  considered  the  possibilit)  that  a  cor¬ 
rection  accounting  for  the  full  finite  electron  beam  width,  and  incorporated  in  the 
starting  e  iirreni  and  elite  ienc\  e  ait  ulaiietns.  e  mile!  e-flee  li\el\  lower  the  theoretical 
second  harmonie  emission  toward  tin  « \perimentall\  incasmed  value’s.  Although 
it  was  felt  that  this  cetrrectmn  was  pi«\s'nall>  more  realistic  than  the  infinitely 
thin  beam  calculations,  ji  die!  tmt  'ucceeei  in  predicting  lower  emission  levels  at 
the  seconel  harmonie. 

The  seconel  consieleraliem  was  that  mode*  suppression  could  account  for  the 
rlerreaseeJ  seconrj  harmonie  moeh  spe-e  trnm  anti  power  levels.  Multimoding  was 


experimentally  verified  between  fundamental  and  second  harmonic  mode*,  a*  well 
ax  between  separate  second  harmonic  modes.  The  quantitative  theory  is  only 
adequately  developed  for  interaction  between  fundamental  modes,  but  qualitative 
theoretical  considerations  support  the  view  that  an  established  fundamental  mode 
will  lend  to  prevent  oscillation  at  2*  .  The  experimental  results  contributed  to 
this  conclusion  t»\  showing  the  strongest  second  harmonn  mode*,  oscillating  in 
gaps  in  the  fundamental  spectrum  It  was  difficult  to  estimate  the  degree  of  mode 
suppression  caused  lo  one  second  harmonic  mode  on  another,  but  the  overmoded 
nature  of  the  cavity  and  ol»*er\ed  2-  2-  mullimoding  definitely  allowed 

the  possihililv  of  this  type  of  competition.  The  qualitative  assessment  based  on 
the  experimental  evidence  was  that  the  fundamental  suppression  of  the  second 
harmonic  modes  was  the  dominant  effect. 

0.2  Rocoiiunendntkms  for  Further  Study 

The  investigation  suggests  there  remain*  a  great  deal  of  experimental  and 
theoretical  atuilv't'  which  could  be  done  concerning  the  harmonic  emission  from 
high  power  gvrofron*  designed  for  the  fundamental  The  experimental  examina* 
'  >*•’,  o’  f  ht  n.f  rroiiM  content  of  'fit  \\  In*  |h  Mlig  gal  let  v  ca\it\  to  compati  vvitli 
rah  illation*  alfcadv  riled  in  this  siucfv  would  provide  an  obvious  first  »lr|i. 

For  the  theoretical  aspects,  the  most  useful  efforts  would  Ice  toward  increasing 
’hr  immediate  applicability  of  ihc  mode  compel  it  ion  theory  to  include  «.•,  *  •  2«t\ 
and  2«»'  •  •  2*  effects.  This  difficult  prospect  would  yield  the  immediate  benefit 
of  quantitatively  testing  the  suppression  which  was  qualitatively  theofired  in  this 
analysis. 
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•ingle  nwdr  an«l\»t*  tiill  alluu*  reasonable  rllinmno  at  n  but  an  rt- 
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